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ABSTRACT 
The directional correlation of the 722 - 556 KEV gamma- 


7 was observed as a function of the physi- 


gamma cascade in Cd 
cal and chemical state of the source. In the case of In€lz in 
a dilute aqueous solution of HCl it was deduced that the undis- 
turbed correlation was obtained, and the following values were 
found for the LeGendre expansion coefficients: A= .090 — ,011, 
A,=.0224 von. It was found that spin assignments of 2-2-0 
for the first three levels of Call* together with the assump- 
tion that the first transition is an Ml, E2 mixture with the 


following percentages: Ml - 95.8, H2 - 4.2 (both + .6), gave a 





good theoretical fit to the observed curve in general agreement 
with results of other workers. The lifetime of the first exci- 
ted level in cail4 was measured and found to be less than or 
equal to 2.3x1071° seconds. To observe the disturbed correlation 
the state of the source was altered to include glycerin solutions 
of several viscosities, the dry salt of InClz, and frozen aqueous 
and glycerin solutions. Calculation of the magnitude of the cor- 
relation perturbing interaction was based on a theory due to 
Abragam and Pound. The liquid sources showed no spoilages within 
the range of viecos ties to which the liquid theory applies: 
however, the dry salt showed an interaction of greater than or 
equal to 700 mcs on the assumption of an electric quadrupole in- 
teraction, and greater than or equal to 165 mcs on the assumption 
of a magnetic interaction in classically describable fields. 

On the basis of field strength estimates in the former case, a 
lower limit to the/quadrupole momentjof the first excited state 


114 ~m24 


of Cd 16 set at .21x10 om’, 
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Ghapter 1 = Introduction 

The first paper on angular correlation of successive nu= 
clear radiations in which the explicit theoretical correlation 
functions were derived was by Hamilton (1H40). It was apparent 
from the begining that here was a very powerful means of inves- 
tigating the low-lying excited states of nucleii. But succes- 
sful experiments had to wait many years for advances in elece- 
tronics to provide coincidence circuits of sufficiently short 
resolving time and for the improvement of radiation detection 
techniques, The first successful attempt to observe an angue 
lar correlation was by Brady and Deutsch in 1947 (1B47)), and 
to Co? 8076 ys | and Cat %* belong the distinction of being 
the first isotopes to have their gamma-gamma correlations mea= 
sured conclusively. But really precise measurements of angue 
lar correlations still had to wait until means of detection 
more efficient and reliable than Geiger counters could be de=- 
veloped. For in this work one measures a statistical function 
and is therefore plagued by statistical errors; hence for pre-= 
cision a great number (hundreds of millions) of quanta must 
be counted. “And under these conditions Geiger counters proved 
quite inadequate. But a new day dawned with the first develop-= 
ment of the scintillation crystal-photomultiplier tube tech- 
nique by Marshall and Coltman (1M47). It was now possible in 
a single day to measure a correlation to the accuracy it would 
have required many years to achieve using Geiger counters, 
And the experimental art was greatly enhanced by the work of 
Hofstadter (1H50) on scintillation crystals of sodium iodide 
activated by small amounts of thallium as an impurity: these 


crystals provided means of measuring gamma ray energies with 


= 





high precision and efficiency. Thus we find in the late 1940's 
and early 1950's a tremendous increase in the work, experimental 
as well as theoretical, in angular correlation, and the first 
paper of Hamilton's is now accompanied by over four hundred 


others. 


are 


Chapter 2 - Theory of Undisturbed Angular Correlation 

The study of the properties of low-lying excited nuclear 
states is greatly facilitated by a nuclear spectroscopic “tool” 
known as “angular correlation", or more completely “angular cor- 
relation of successive nuclear radiations", The short lifetime 
of most excited nuclear states makes their investigation by die 
rect means such as resonance methods next to impossible. Hence 
we must probe these states indirectly, and the method that comes 
immediately to mind in the case of low-lying states is to study 
the radiations which preceed and succeed them. When we investie 
gate some property of these radiations as a function of the angle 
between their propagation vectors we do "angular correlation", 
As one can imagine, the field is very broad; however, it can be 
subdivided both as to the nature of the radiations involved and 
as to the particular property of the radiations under investiga~ 
tion. In the former case, the theory is broad enough to ine 
clude any of the well known nuclear emissions, viz. beta rays, 
gamma rays, alpha particles, neutrinos, in any allowed combina- 
tion, vize B-¥,\X-X,<-¥,A-@,@-v. In the latter case, there 
are two general properties of the radiations which divide the 
field of angular correlation into (a)) "directional correlation", 
wherein one is interested in the relative probability of emise 
sion of the second radiation as a function of angle relative 
to the first, and (b)) "polarization correlation", wherein one 
is interested in the absolute or relative polarizations of 
either or both of the radiations as a function of included angle. 
The original theoretical investigation in this field was cone 
cerned with gamma-gamma directional correlation and as mentioned 
previously was done by Hamilton (1H40). Since this investigation 


is concerned with gamma-gamma directional correlation, we will 


an 





sketch Hamilton's arguments in some detail. 

Hamilton starts by considering a nucleus in an excited 
state, Ao , where 2 refers to the magnetic quantum number and 
A refers to all other quantum numbers. This nucleus is to 
decay to a ground state, CG , by way of an intermediate state, 
B . In so doing it must emit two gamma quanta, each characte= 


rized by a propagation vector, K » and a polarization unit vec 


tor, eS ; (The propagation vector, K= K Ko » where Ke is @ 
unit vector in the direction of propagation and K= oe where 
\ is the wavelength of the gamma ray). Subscripts | and 2 


on the propagation vectors hereafter will refer to the first and 
second transitions respectively. 

There are two very general propositions concerning a radia- 
ting system of such nmucleii which underly this theory. They 
are stated below without proof: 
(1)’ Assuming the radiation is detected by direction insensitive 
counters, each of the three states involved is oriented complete- 
ly at random. The expressiom “oriented completely at random" 
as used here has a very precise meaning, which is: 

(a) equal populations among the magnetic substates, and 

(bd) random, i.e. unrelated, phases among the wave functions 

of the magnetic substates, 

(Orientation at random is necessary in order that each substate 
be an independent radiating system)’ 
(2) Each of the two readiations involved is completely isotropic. 

At first it may seem nothing short of amazing that out of 
propositions (1) and (2), which apparently express the "acme of 
isotropy”, we can derive a theory which predicts an anisotropy, 
iJe. a correlation. 


a Wa 





But this is so, and this idea occupies a central place in the 
theory. The clue to resolving the apparent anomaly is that if 
we put a condition on the population we consider and count only 
that sub-population in B whose preceeding radiation went in @ 
certain direction, 0, , with respect to the z axis, then the dis- 
tribution of this sub-population is not necessarily equal nor 
are its phases necessarily random This condition can be come 
pared to that reported by the census taker, "The population of 
this state is uniformly distributed in area, however those whose 
parents came from Ireland live predominantly in the northeast", 
So much for similes. 

This idea is important enough to belabor a little further 
with a quantitative approach. Let the relative transition pro- 
bability per unit time from state Ax to state im be ems 
Then from (1) the total rate of radiation from each substate in 
A must be equal, hence: 

(3) a Q aon is independent of & , and since the various sub-= 
states of BK must be equally populated: 

(4) 2 4am is independent of M . And since the radiation 
is isotropic by (2), each mode of a given multipole order of ra- 
diation must be equally excited:: 

(5)) 2 99,000 is independent of mm , where —L <mm<<+L if L 
is the multipole order of the radiation, 

Let Pr,mlo\be the angular distribution of a multipole of or- 
der L carrying an angular momentum whose z axis projection is 
given by am, The relative number of quanta in the As—->Ba. trane 
sition for example that proceed into a solid angle du at an 
angle of 6, with the z axis is given by am $1, 2-m (B,))dus, ry 
ean further be shown that the following is a property of the 
distribution functions? 





(6) x: dim(O.) is independent of & . 

We will now proceed to show that these relations imply ma- 
thematically that the first transition is isotropic ~ - - which 
they must, for they were based on that hypothesis (2). 

From the conservation of angular momenta: 

im, 2 =n + om 

Let the angular distribution of the first transition be 
given by #\(6,) . Then: 

(2) Bar = ZF 2 Ae.d-em Pr, wn (On) 

Applying (5), (6) becomes: 

(9) @cay = 2 ae,e 2 Prism lO 

which by (6) 46 independent of ‘oF » nence isotropic. 

We will now show mathematically that the total population 
ef the states Nm is independent of ™ , which again is directe 
ly from hypothesis, but that the sub-population, arising from 
radiations in the direction 6, is not necessarily independent 
of m , hence the second radiation from these particular sub- 
states is not necessarily isotropic:: 

Let Pm be the relative probability that a state Dm is 
occupied. . 

Let Pon (01) be the relative probability that a estate Ben 
is occupied through a radiative transition proceeding in the 
direction 6; with respect to the z axis. 

Then we can write immediately: 

(10) Tm = ad Qiao which by (4) is independent of m , and:: 
(11) Pm(6r = 2 Qos gem (Bilwhi.ch by (10) is not necessarily 
independent of qm . This was to be shown. 

After going through a derivation based upon a system of 
quantized radiators and field oscillators, Hamilton obtains the 


following expression for the complete correlation function: 


> Bs. 





ee oe * | 2 
(12) W (Kei. Kor) = = = Wa 2 (Nel Hi Kor etl Bl Bm] HRon8a)| CP) | 


ce 2,P 





where Cy and - represent averaging over all non=observed 
f 
properties of the radiations | and Cd such as polarization, and 
where H (Ho,e » is the perturbing Hamiltonian, given by: 
li =| Ke 
* = Ae —Gaec C- e, dv 

(12a) Hike) = - \i Adv % \e 

where A = vector potential of the radiation field taken 


in the gauge where scalar potential is zero, 


= 
L 


current density of the radiating system, 


charge density of the radiating system, 


a ao 
TT 


time dependent amplitude of radiation field oscil- 
lator and is independent of position, 

n= position vector, and 

© Ko, as before. 

The cross terms in (12) resulting from the summation over m 
before squaring make it very unwieldy. Hamilton goes on to show 
that these cross terms arise because the phases in state B are 
no longer random, that if the phases were random the cross terms 
would vanish, and hence the summation over ™ could take place 
after sqnarink. He further shows that the phases are random if 
one takes the z axis as the direction of the first emission. 
Several authors have proved this theorem also, but an especially 
simple proof is due to Lippmann (1L52) and is included here be- 
cause of its physical insight. Lippmann's proof is based on a 
fundamental quantum mechanical theorem: “ — the probability of 
a transition from an initial to a final state, via a set of in- 
termediate states contains interference terms only when it is 
not possible to measure ( — without disturbance - ) the inter- 
mediate state; if the experiment is inherently capable of spe- 


cifying the particular intermediate state the system passes 
= 2? 





through, the classical law of compounding independent probabili- 
ties holds and no interference terms occur -", Since we are 
measuring precise angles, a plane wave description of the emitted 
particles is appropriate ( - re and Rene well defined ~ ). Now 
if the z axis is chosen parallel to es, » then the description 
of the emission of the first particle is invariant under a rota- 
tion of the coordinate system about ry (assuming we average 
over polarizations): hence the total magnetic quantum number as 
well as the quantum numbers of the emitted particles (Liam, ) 
are constants of the motion. (It can be shown by other means 
that under these conditions ~m,=2 | )< If we consider our ra- 
diating system to be a single nucleus, the substate, Aw » before 
the first emission is in principle measurable, and since our mea~ 
surement of the first quanta did not disturb the system, the 
specification of substate sige follows on the basis of the con~ 
servation of angular momentum. Then the theorem (about the ab- 
sence of cross terms) follows. 


Thus when we define ® as a unit vector parallel to the 


z axis, we lose our cross terms and (12) reduces to? * 
a3) W@ jk.) = a sas, ts (a 1H&29IBe) Belfi gebcdl 
8B) 


If we let the relative probability for a transition from 
state Rm to state P ata direction 8 relative to the z axis 


be Pre (8), then: re 
(14) Pw (oe) => maine | (Ce \Hckor ea Bm) 
Z 





And with a similar aefinition for the first transition (13) 


becomes: 


as) W(e) = te. Pam (0) Pre (6) 


Ry Mm ,F 


@- 8 = 





(15) is presented here because of the interesting manner 


in which the individual transition probabilities decompose: 
~ VY) 2 
CRI camera is len lea cians lili sia lara ca lie 


with La Wee 
and Mm = Mat P 


where [bh = the spin of state Rm, 
lo = the spin of state Coe ’ 
L2t= the multipole order of the second radiation, 
MM = the magnetic quantum number of state Bm, 
*¥ = the magnetic quantum number of state Cp » and 
WN. the magnetic quantum number of the second 
radiation. 

The first factor in (16) is the well known Clebsch-Gordon 
coefficient in the notation of Condon and Shortley (1635). It 
depends only upon the quantum numbers of the radiation and the 
states involved, and gives the intensity for the particular 
transition. The second factor gives the angular distribution 
for the particular transition and depends only upon the quantum 
numbers of the radiation. Falkoff and Uhlenbeck (1F50) give 
the general expression for the angular distribution functions 
and the particular examples given are for diple and quadrupole 
radiation, 

Hamilton expressed (13) finally as a cosine power series 
in which due to the properties of the Pain only even powers 


appear: 


K! 
(17) WW (Q) = } Oomcos** 9 
K-@ 
The series terminates and Yang (1Y48) showed that the highest 


K is given by: 
(ja) K’SLI, K’S La | KS Tt 


) 


whichever is least. 





Hamilton calculated °%%/o, and “%‘4/q, for all combinations 
of multipolarities up to and including quadrupole. 

It was later found that if one expresses the correlation 
function as a tories, in LeGendre polynomials, 

(is) Wie, = = Nx Pex (case) 

the coefficients, i can be expressed quite generally and 
explicitly, and that (18) still applies. This was first shown 
by Gardner (1649) for the case of conversion electrons, and was 
later extended by others to arbitrary radiations. It turns out 
that the Azk can be decomposed into two factors: 

(200) Azxw = CARY2K (RC) 2k 

where the first factor contains parameters relevent only 
to the first ( A—~-B ) transition and the second factor, those 
parameters relevent only to the second (K—-C ) transition. 
The expression for the first factor is given below, the second 
being similar in all respects: 

(021) (AR)ax = ZL Limmi-omil Labs KIC) (en 
OWN GT ae De 

where the W is the “Racah coefficient", Racah (1R42), and 
all other factors are as given before (16). 

Using (21) we can now calculate gamma-gamma directional 
correlations explicitly, with a few reservations as indicated 
below. For example, one can find the Clebsh-Gordon coefficients 
in Condon and Shortley (1035), the Racah coefficients can be 
found in (2B52), and the ee (0) can be found in (1F50)<: 

However, these correlation coefficients, Az« » have 
been calculated explicitly and presented in convenient form 
by Lloyd (1L51). 


Now, it seems, it would be appropriate to point out two 


=] 10 «~- 





limitations on the theory as developed so far and as applied to 
gamma=-gamma directional correlations: 

Limitation 1 —- The above theory refers only to pure multi- 
pole transitions. It cannot handle, for example, parity forbid- 
den transitions which proceed partly by magnetic 2"_pole and 
£1 


partly by electric 2” -pole radiations. 

Limitation 2 = The above theory assumes that the interme- 
diate state is unperturbed. As was pointed out previously, the 
existence of a directional correlation depends upon unequal popu- 
lations among the various sub-states Km. A perturbation ap- 
plied to the intermediate state would tend to cause a mixing 
of these states and hence a relaxation towards a uniform popu- 
lation. Any such equalization would then by propositions (1) 
and (2) above cause a reduction in the correlation. 

Limitation 2 has been removed by Ling and Falkoff (1L49). 
They calculated the effects of various mixing ratios in the mag- 
netic n-pole and electric n/l-pole mixed transitions. Their 
results however are complicated by the presence of an undeter— 
mined relative phase factor in the ratio of these matrix elements. 
Lloyd (2L51) showed however that only a relative phase of O or 
180° had physical significance. 

Limitation 2 is still very much with us. However it has 
turned out to take on more the qualities of an asset than a 
limitation. For by studying the reduction in the correlation 
caused by the perturbation of the intermediate state we can 
learn something about the electric and magnetic moments of 
this state. This study becomes the field of “disturbed angue 


lar correlations" and is the subject of the next chapter. 
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Chapter 3 = Theory of Disturbed Angular Correlation 

General interest in the subject of disturbed angular cor 

relation arose when several correlations were observed which did 
not fit any consistent theoretical predictions. The situation 
was further confused in that the observed correlations seemed 
to depend on the physical and chemical state of the source. 
The most famous of these “anomalous correlations" is the gamma- 
gamma cascade in tae the daughter by K capture of . ae see 
for example Boehm and Walter (1B49)), Roberts and Steffen (1R51), 
and Frauenfelder (2F51). 

It became apparent that these "“spoilages" of the correlation 
might be the result of a perturbation of the intermediate state 
during its period of existence between the emissions of the 
first and second quanta. For, as we have seen one of the basic 
assumptions of the “undisturbed” theory is that the intermediate 
state remain unperturbed during its lifetime. Several authors, 
including Goertzel (1G46) and Alder (1A52) have developed suc- 
cessful theories based on the interaction of various nuclear 
moments of the intermediate state with perturbing fields. We 
shall confine ourselves to tracing a theory due to Abragam and 
Pound (1A53) which is recommended by its generality. 

Abragam and Pound show that the effects of such a pertur- 
bation can always be represented in the form: 

a) Wie) = a (S 2k Ak (aK (coee) 

where the Azkx are the “unperturbed” coefficients and the 
(Szx which contain all of the perturbing effects act as damping 
factors on the unperturbed coefficients. 

They point out several interesting points in connection 


with the (Fox relative to the type of source; namely: 


Sie = 





(a) If the source is in the form of a polycrystalline 
powder the Gex approach a minimum or “hard core" value as the 
magnitude of the perturbing interaction increases without limit; 
hence the correlation is never completely destroyed. 

(b) If the source is in the form of a single crystal, 
O< (5.x <.| » Gepending on the orientation of the crystal axis 
and the strength of the interaction. 

(c) If the source is in the liquid state the (>2x can 
approach zero as the magnitude of the perturbing interaction 


increases without limit, 


Abragam and Pound start with Hamilton's basic directional 
correlation equation (2-13), slightly modified to conform to 


the former's notation: 


(2) Wie) = S 2. CalnleolHele\clHelG/ WGI la) 


1,2 Ay,C, jC’ 
where O,C represent quantum numbers of the first and 


third states respectively, and BAe represent quantum numbers 
of any complete set in the intermediate state, Id ° H, ana H2 
are the interaction Hamiltonians for the emission of the first 
and second particles respectively (See 1-12). 

Let us suppose that the first emission takes place at 
time Lb =O, and that during time t the intermediate state is 
perturbed by some Hamiltonian “K « The effects of this inte- 
raction can be represented by an evolution operator, LI) (+) 
acting on the intermediate state: 
i <(t) = Uc 

where X(t) represents the intermediate state vector after 
transormation by the interaction. 


——— 


 ) Ce) can be represented as follows: 


» OS = 








a 4 
Gay Wa) — e h 
ase. S 
for static K , or: 
pam de (Edt! 


(4p) Ua) = © 


for time dependent K , 
If now we substitute (3) into (2) we have the time depen- 


dent correlation functions: 


(5) avty=GB2_ (al | QI D%erHale cle Ueel6G 1H. la) 
where Ut he ap ae re conjugate of TT ; 

We must now consider the probability of emission of the 
second particle between time t and tidt., This was ignored 
in the undisturbed theory because it yielded only a constant 
multiplicative factor. However, in our case we must consider 
the fact that the relative probability for emission of the 
second radiation between time t and t+dt is given by: 
lege en 45 
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where we consider the source as a whole. If we record all events 





taking place until time U2 , our function becomes: 
to 
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Ir K is static we can choose for GB , the eigen vectors 


of K , which we define as 1b ; then (6) becomes: 
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and if we let t2—-co , we have finally the result ob- 


tained also by Goertzel (1G46):: 
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Abragam and Pound's next step is to reformulate the corre- 
lation function in a form convenient for the isolation of the 
perturbing terms. They use an expression similar to one derived 


originally by Alder (1A52): 
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where the VER) are ordinary spherical harmonics in 
which [L, and (2, specify the direction of emission of the first 
and second quanta respectively, 1. and il are coefficients 
which are independent of the perturbing Hamiltonian; and ali of 
the effects of the perturbation are contained in eee: 

If the perturbing field can be described classically, 


Uke wut) is expressed by: 
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where the first two factors are Clebsch-Gordon coefficients 
in the notation of Condon and Shortley (1035), the U is an evo= 
lution operator, and the m’s are the eigen values of the inter- 


mediate states Ban which are chosen to be eigen functions of Li. 


If the perturbing Hamiltonian, K , is static and if we 
define 1 as an eigen state of this Hamiltonian, then the U 
can be represented as: 

‘ 1 iat 
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With this substitution (10) becomes: 
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Now, considering the source as a whole just as in (6) if all coin- 
cidences are registered we must integrate U’(t) with a decay weigh- 


ting factor: 
aad — 7 es Pra. 
i ‘TN 
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Then the integrated coefficient, I(qupipids becomes: 
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From (9) and (13) we can see why . the interaction Hamil- 
tonian for the perturbation is invariant under a rotation about 
the z axis and if the first particle goes in the direction of 
the z axis then the correlation is undiminished. In the first 
place, if the particle goes in the z direction, the only term 
of the Yn. (0) which does not vanish will be where =O 
Since the Hamiltonian is invariant under a rotation about the 
z axis, it commutes with Ll< » hence Lb= Emand from the first 
Clebsch-Gordon coefficient in (13) we see that since (™M =im’‘+ 41. 
that omen’ «. Thus Eb=Eh! in (13) and the denominator re- 
duces to unity. Now due to the various summations over lb ané 
Lh’ » we have the multiplicative factors of Qem orn" and dem” con! 
which when we remember that an=cm’ , imply that m= am” 
hence /{t;=Q. Then nine Ck, Kk, OG) = oan ie Ke 
wnich ig the normalized form. We have thus shown that under 
these conditions the theory predicts the well known experimental 
result that the correlation is unperturbed. 

We are now ready to consider the problem of the crystalline 


powder, or polycrystalline source, as it is called. The observed 


correlation, W(¥), ig obtained by averaging (9) or (13) over 
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all possible orientations of the various micro-crystalline frames 
of reference which is equivalent to integrating the spherical har- 
monics, UCN pate » over the unit sphere, keeping always the 
angle between Ni and str fixed; this angle we call Y . If we 
do this we find from the orthogonality of the spherical harmonics 
that the only non-vanishing products are where ki. =Ki= kK = and 
H.=p_r. = fe ,. It also turns out that the integral so defined is 
independent of (L . The observed correlation can then be writ- 


ten, if we average (9) to obtain the time dependent correlation: 
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If we average the integrated correlations we obtain? 
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where 
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Notice that the (5x factors which contain all of the effect 
of the perturbing Hamiltonian are simple reduction factors which 
cause no mixing of the unperturbed coefficients. Also, the CG 
factors depend only upon the intermediate state and not upon 
the preceeding or succeeding radiative transitions. 

We give below an explicit expression for the time depen- 
dent as well as integrated Gx'swhere the perturbing field can 
be described classically and has axial symmetry with respect 


to the micro-crystalline axes: 
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So far we have said nothing of the nature of the perturbing 


interaction, ate; whether magnetic or electric in nature. If we 
assume for a moment that we are concerned with the interaction 

of the nuclear electric quadrupole moment with the average elect- 
ric field gradient, classically describable and axially symmetric, 
at the nucleus, then the perturbing Hamiltonian is already diago- 


nalized and its 0 ey values are; 

CS a me ay 

4T(2T-1) 

where CQ is the nuclear quadrupole moment, defined: 


eQ = i; \e (Scos* Os-Nny VY sa At 
321 
the summation going over all protons where (Qj is the angle 
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between: and the nuclear symmetry axis. 

If we now substitute (20) into (18) we obtain an explicit 
expression for the effects of the electric quadrupole interaction 
which for I= 2 we write: 
in) Cy = (= + ae eas & Ecos suet + 2 C054 Wee t | 
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where Uk, is the energy difference divided by / between 
the levels which are split the hs the interaction; for 
me WS az 
I= 2 this yields, “oe. = — 

If on the other hand we consider a perturbation due to a 
classically describable, axially symmetric with the various 
micro-crystalline axes, magnetic field which interacts with 
the magnetic moment of the nucleus, we have: 


@ Ca =—e- kh = - Heap = - hwo, Hom 


where {Jt is the nuclear magnetic moment in nuclear mag- 


ay A 





netic moment in nuclear magnetons, H is the magnetic field, am 
the magnetic quantum number, and +», is the angular velocity of 
Larmour precession given by pH ° 


AL 
As before, b3 introducing (22) into (18) we obtain for I= 2: 


(23) ena ==(! 4+ 2cosus t + @cos % Wem O | 
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There are two important features of (21) and (23) which 
bear noting: 

(a) The G, (it’s are periodic functions of time with a maxi- 
mum value of unity and a period of — - This periodicity arises 
from the axially symmetric nature of the fields which split the 
energy levels in rational multiples of the minimum splitting. 
Fields of lower than axial symmetry remove the necessarily pe~ 
riodic nature of the G,it)’s. 

(b>) The integrated Gx (ts)’s , given by (17), which are 
observed by a coincidence circuit whose resolving time is long 
compared to the nuclear lifetime, have a lower limit (given by 
the leading terms in (21) and (23)) known as the “hard core", 
This implies that the anisotropy will also have a “hard core” 
no matter how strong the interaction, wets . As will be seen, 
this is not necessarily so in liquid sources, where a sufficient- 
ly strong interaction can reduce the anisotropy to zero. 

The theory of spoilage by axially symmetric, classically 
describable perturbing fields in polycrystalline sources as 
described above has been tested experimentally and in some cases 
(1K53) gives correct results. However, the results of Aeppli 
et. al. (1A51) and Steffen (2853) (1853) on various compounds 


ter ALLA 
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of Cd the daughter by K capture of In show reductions 


in anisotropies to values lower than predicted by the theory of 
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axially symmetric perturbing fields as put forth above. The un- 


disturbed anisotropy of the gqiil 


cascade is about -0.2; the 
axially symmetric theory predicts a “hard core" anisotropy of 
about -0.05. For example Steffen (2853) reported among others 


the following anisotropies: 


Source:: Anisotropy: 
In€lz dry @ 20°C -.012 + .005 
@ 540°C -.022 + ,.006 
InI, dry @ 20°C -.020 + .006 
@200°° -.021 + .006 
Average of above -.019 + ~ 003 


It is not surprising that the spoilage theory based on the 
assumption of simple axially symmetric fields does not give 
precise results. We shall present below a discussion of two - 
of the possible mechanisms for the spoilage of polycrystalline 
sources below the predicted hard core values. 

The first explanation has to do with the symmetry of the 
fields experienced by the nucleus. The fields in polycrystal- 
line sources are extremely complex, and even if the indium a- 
tom in a single isolated crystal were in a field of axial syne 
metry, when there are many millions of microcrystals close to- 
gether, oriented at random, and pressing on one another from 
all sides, there is good reason to believe that the aforemen- 
tioned axially symmetric fields could be distorted to a lower 
degree of symmetry. Also, following a process such as K cap- 


ture in tnitt 


» the daughter cadmium atom may have sufficient 
recoil energy to displace it from its former axially symmetric 
lattice site. Abragam and Pound (1A53) calculate the effects 


of fields of lower than axial symmetry. They give the hard 
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core values for Gx’s in rhombic fields for various values of 
the spin of the intermediate state. They find that there is in- 
deed a reduction of the Gx’s over and above that obtained with 
axially symmetric fields; however, the resulting reduction in 
the anisotropy is still not sufficient to explain the results 
of 1853 and 1A51. For information we list here the results of 
their calculations for the hard cores for I= 2: 
Rhombic Fields Axially Symmetric Fields 
Go(1im) 10/35 13/35 
G,(1im) 18/63 29/63 
It may be of interest to note here that if the pertur- 
bations remove all degeneracies, the hard core value of G. 


should be G,(lim) = Comparing this with the above 


results we see that the rhombic fields apparently remove some 
but not all of the degeneracies present in the axially symme- 
tric case. 

The second possible mechanism for spoilage below the pre-= 
dicted hard core value, and it seems the most fruitful to date, 
has to do with the excited state of the electronic configuration 
of the daughter atom. Immediately following K capture the daugh- 
ter atom is in a very highly excited atomic state with a hole 
in the innermost shell. The electronic configuration returns 
to the ground state by the emission of X-radiation and the pro- 
Guction of Auger elections as the “hole” moves outward. As 
more and more Auger electrons are produced the atom tends to 
assume a considerable charge, and each Auger electron itself 
adds another “hole” to the configuration. After a short while 
the atom begins to “look” like a "Swiss cheese" as Frauenfelder 


puts it. Under such conditions the nucleus experiences extremely 
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large and rapidly fluctuating magnetic fields and electric field 
gradients due to its own electron cloud. It seems possible that 
this could cause a time dependent perturbation very similar to 
that which would be experienced by the nucleus of an ion ina 
liquid under the influence of the random motion of its neighbors; 
let us call this the “confused configuratiom effect". And in 
a liquid, as will be seen, there is no lower limit or hard core 
to the anisotropy for a sufficiently strong interaction. On the 
other hand if the nucleus were to "see" during its lifetime 
fields due to its own electrons which were constant for a while 
and which would then vary stepwise, it is possible that it could 
get down to one hard core value in this time and then with the 
stepwise variation in the electronic configuration, start down 
toward a new and still lower hard core. Let us call this the 
"quasi-static configuration effect", 

To test the validity of the K capture excitation hypothesis 
Kraushaar and Pound (1K53) performed a directional correlation 


experiment using the 396 KEV isomeric level of gqitt 


as the pa- 
rent of the decay. Here, the electron shells of the cadmium 
atom are unperturbed during the decay, except possibly by ine 
ternal conversion. Kraushaar and Pound used several sources 
differing in chemical and physical state. It was found that 
the maximum spoilage of the anisotropy occured in a polycrystal- 
line source of hexagonal CdClo, and significantly this reduc- 
tion was down to, but not below, the hard core as predicted by 
the theory of axially symmetric fields. 

We will now introduce what is called the "shell relaxa- 


tion time", C, , which we can think of as the mean time required 


for the electronic configuration to return to the ground state 
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following a strong excitation such as K capture. We list below 


two typical values of 2s for different materials from (1H52): 


Material: : 
Impurity Centers in yy 107 "seconds 
Ionic Crystals 
Metals ~ 107) “seconds 


We can fit three experiments into the pattern of explana- 
tion of the “excited configuration" hypothesis. In the experi- 
ment of Albers-Schonberg et. al. (2A53), using a source consis- 
ting of a single axially symmetric crystal of indium metal, the 
effect is noteworthy because of its absence. For certain orlenta- 
tions of the crystal axis the observed anisotropy was essentially 
unperturbed which hardly seems possible on the basis of the “ex- 
cited configuration" hypothesis. This can be explained, however, 
as the result of an extremely short shell relaxation time in 
metals (~ 1071 seconds, above). In cqlil the lifetime of the 
intermediate state is 8.5x1078 seconds, and when the shell re- 
laxation time is very much shorter than the nuclear lifetime the 
effect has too little time to take place and shows but little 
disturbance to the anisotropy. This is similar to the case in 
liquids wherein the correlation time is very much shorter than 
the nuclear lifetime; the liquid theory will be presented later. 

On the other hand, the experiment of Steffen (2853) cannot 
be explained om the basis of the “confused configuration" effect, 
if we believe the figure’ above of 2, Sy 10 8s. for impurities in 
ionic crystals. For here 73; >tn, and under such conditions the 
nucleus sees practically constant fields due to its own electro- 
nic configuration, as will be seen. This requires the application 
of a modified polycrystalline theory considering the quantum 


mechanical coupling between the nucleus and electron cloud, taking 
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into consideration the fact that the angular momentum of the 
shell, J, is changing during the nuclear lifetime. This, it 
would seem, is an example of the “quasi-static configuration 
effect" mentioned earlier. It may be that the reductions in 
anisotropy below the hard core values observed by Steffen (2853) 
could be due to a combination of two independent effects of com- 
parable magnitude, viz. (a) the "quasi-static configuration ef- 
fect", and (b) the crystalline field effect treated by Abragam 
and Pound (1A53). In a situation like this we could have an e- 
lectric quadrupole interaction with the crystalline electric 
field gradients and simultaneously the superimposed magnetic 
and electric interactions with the "own-electron” configuration 
which have not yet "cooled" down to the ground state. The 
theory has not yet been worked out for such a situation: but 
this seems a reasonable explanation for the occurence of Gy 
hard cores lower by a factor of two than those predicted on the 
basis of "axially symmetric fields" alone. 

Futher justification of the "quasi-static" view follows 
from a simple estimate of the times involved for the outward 
progression of the K-capture produced hole in the atom. For 
a first approximation, if we consider the atom alone independent 
of its neighbors and neglect the additional holes produced by 
Auger electrons, then we can calculate the time required for the 
radiation process on the basis of a model in which an electron 
jumps from the n shell to the n-1 shell with a speed given by 
shielded hydrogen wave-functions. The first transitions between 
highly excited atomic states of large energy differences are 


9719 


very fast, being of the order of l seconds; while the later 


transitions are relatively slow, lying in the range for optical 
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=8 seconds. But the important point is 


transitions, or about 10 
that the times for the various transitions are related roughly 
as a geometric progression, the second requiring about 10 times 
the time required for the first, etc. This means that on the 
average during the lifetime of a particular nucleus it has ex- 
perienced essentially a constant field resulting from the last 
electron configuration in existence before the nucleus decayed 
to the final state. In addition, the sum of the times required 
for the early transitions , if less than on seconds would 
cause no perturbation just on the basis of its short duration. 
All of the above seems to lend support to the “quasi-static hy- 
pothesis as the additional spoilage mechanism, 

On the other hand, because of the possible similarity of 
the processes experienced by the nucleus after K capture to those 
in a liquid, we will sketch the results of Abragam and Pound's 
(1A53) theory of correlation spoilage in liquid sources. 

The time dependent correlation in liquids is analyzed using 
(9). Taking the z axis as the direction to the first counter, 


fei vanishes as discussed under (13), and since there is no 


preferred direction in the source Jtz also vanishes. Then we 


have: 
° wo 
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where [(no‘\Uitlm) | is the probability that a nucleus in 

an initial state m at t = O will find itself in state m’ at 

time t. This relaxation process is what spoils the correlation, 


for as pointed out in Chapter 2, an unequal population among 
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the various magnetic substates of the intermediate level is ne- 
cessary to produce a correlation. Furthermore it can be shown 
that the rate at which this relaxation process takes place is 
related linearly to the rate at which the correlation coefficients, 
Ao,, (1) diminish. 

In the development of the theory the concept of “correla- 
tion time", CT. , enters. Qualitatively, one can think of the 
correlation time as the mean time required for a nucleus to in- 
terchange its neighbors. 

It turns out in the development of the expression for the 
time dependent transition probability above that when f?7tc, but 
with the condition that the transition probability is very much 
less than unity, we can use the theory of polycrystalline sources 
to describe liquid phenomena. We apply this result by expanding 
(21) in a power series in t and replacing tb py 2tte . 

It is interesting to note that although the G,'s thus obtained 
fall off as the first power in t, they are still larger than 
the G's from polycrystalline sources because tf-Kt . Similar- 
ly the not unexpected result that where T.-7t we may apply 
polycrystalline. solid results to liquid sources evolves from 
the theory. This is reasonable in that the fields of an in- 
finitely viscous liquid would approach roughly those of a solid 
except perhaps as related to spatial regularities. It is also 
interesting to note that in this latter region the "fall-off" 
of the G,'s as a function of time is parabolic as t—> 0, just 
as we find in (21). 

In the region in between where the correlation time is 
neither large nor small with respect to t, explicit results 


are developed applying only to liquid sources. These results 
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apply also where @<t , but not where 7 >7tT . The integra- 
ted G's for I= 2 are given as follows: 


Electric Quadrupole Interaction: 


| | 
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Magnetic Dipole Interaction: | 
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It may be of interest to compare the shape of the spoiled 
anisotropy curve based on G,'s which come from no particular 
physical model with that resulting from (21) and (23). The 
purpose is to show the similarity in shape for small values 
of Wotw of the anisotropy curves arising from Gy(t)'s quite 
Gifferent in nasure. This is done in Section 5.4.4 based on 
the following remarks: 

Let us construct a &,(t) in the form of a periodic square 


i = = 
i 
i 
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wave of half width %« » as shown 





in the diagram at right. Then the 





integrated G,(wetw)'gs which we will 





label Gy} (Wotw ) are given by: " ¢ josh == 
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Note that the hard core values are given by bk 
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The Gy can be related by applying a modified form of Abragam 


and Pound's Eq. 41 (1A53), for integral I: 


(29) 7 (4k4l) dex = (2040 7 


K=0 


EO Oia 





Chapter 4 ~- Choosing an Isotope 
Out of all the hundreds of isotopes, how does one pick one 
on which to perform an angular correlation experiment? We sketch 
below the procedure we used in the hope that it will be of some 
future assistance. The procedure consists in giving each "can- 
didate" a series of tests and using the comparative overall re~ 
sults as a basis for selection. 
These tests can be classified as follows: 
1 Decay Scheme 
2. Availability 
5. History 
Decay Scheme 
The first thing that must be done is to establish the exis- 
tence of a cascade. Then the disturbing influence of radiations 
within the energy resolving power of the apparatus must be con- 
sidered, and whether these radiations are also in cascade. It 
will be found in connection with “History” that most of the simple 
cases have already been studied, hence the necessity for good 
energy resolution will become more and more important. If we 
consider gamma—gemme cascades it will be found that in connec- 
tion with the energies of the radiations there are two very im~ 
portant considerations, viz. (a) that if the lower energy is 
above a certain minimum then one can use “energy selection shiel- 
ding" as discussed in Section 5.2.1, and (b) if both energies 
of the cascade are close enough to one another so that each 
channel can be sensitive to both radiations, then the coinci- 
dence counting rate will be up by a factor of four over the 
case in which each channel is sensitive to only one of the 
radiations. (1H53) is a good reference for decay schemes and 


various other information concerning isotopes. Interpretation 
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of results becomes easier when we know that at least one of the 
transitions involved is a pure multipole as we saw in Chapter 
2; this is insured if we choose an even-even nucleus where in 
all known cases except one (Ge’“) the first excited state has 
a spin-parity assignment of 2=-plus and a ground state of zero- 
plus; this insures that the transition between these states will 
be pure electric quadrupole. See (3853) for information on 
excited states of even-even nucleii. The total internal con- 
version coefficient is also a limiting constraint in the choice 
if one wishes to do gamma-gamma angular correlation. We chose 
an arbitrary upper limit to the K-plus-L coefficients of 10: 
anything above this number we would reject. In general it turns 
out that for a @ greater than 60 and energy less than 175 KEV 
the transition will be too highly converted for E2 transitions. 
(3B52) page 618 and (1G51) give formulae and experimental re- 
sults relating to internal conversion. As was seen in Chapter 
3, in an investigation of “spoiled” correlations the lifetime 
of the intermediate state is important, and in general should 
be greater than ome seconds, although this lower limit is 
somewhat flexible. (1451) contains formulae and experimental 
results relating to lifetime. 
Availability 

The availability of an isotope is determined by the availa- 
bility and lifetime of the parent. If local facilities are unable 
to produce the isotope then we are limited in general to parent 
lifetimes of at least several days. (1H53)) has information re- 
lating to the production of isotopes. 
History 


Unfortunately there is as yet no comprehensive survey of 
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experimental results in angular correlations under one cover. So 
after passing "Decay Scheme" and "Availability" tests and isotope's 
angular correlation history must be looked up. Of particular in- 


terest in disturbed angular correlation work are those cases 


where previous work has shown anomalies, such as in par, 


gil4 


Our choice of ¢ was based on considerations similar to 


the above and also due to the history of spoilage in its sister 


isotope eqttt, See Fig. 4-1 for the decay scheme. Our interest 
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lay in the first excited state in Cd via the 722-556 KEV 


gamma-jfamma cascade. 
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Chapter 5 = The Experiment Part 5.1 — Apparatus 
Section 5.1.1 — Apparatus, General 

& directional correlation apparatus must provide means for 
the detection and exclusion of radiation and for the recording 
of events. The term "directional" brings to mind the question, 
"Relative to what?" The answer necessitates the use of two 
channels, one to fix a direction in space, the other to measure 
a relative angle. 

The exclusion function of the apparatus is very broad. 

First of all it must exclude radiations proceeding at other 

than the selected angle. Second, it must exclude in each chan- 
nel pulses arising from unwanted radiations of various energies. 
Third, it must compare the pulses in the two channels in time 

and exclude all those which do not occur simultaneously or within 
a@ small time,Z , of one another, 

Finally, the apparatus must record certain unexcluded, or 
selected if you will, events occuring during a given time inter- 
val and present them in a convenient manner. To accomplish these 
ends the appardtus is divided into mechanical and electronic 


ra 


components. < 


Section 5.1.2 — Apparatus, Mechanical 

The mechanical parts of the apparatus are extremely simple 
and consist of the source box and its included angle measuring 
assembly. 

The source box provides means to mount the source and angle 
measuring assembly, absorb unwanted radiation, as well as shield 
the immediate neighborhood from the soft components. It also 
provides means to change or maintain the temperature of the source 
and to inhibit frost formation on the source holder, which is 


accomplished by keeping the box airtight and drying the entrained 
«= 62 . 
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air. (Activated Alumina 8-14 mesh was used as a descicant and 
performed satisfactorily). The source box measured 33x33x16 
inches; it was designed with a removable top with an additio- 
nal opening for easy access, normally closed with a lucite 
window. A small hole was provided in the center of the top 
to serve as a liquid air filling connection and another hole 
was provided in the lucite window for dry ice filling. 

The angle measuring assembly (Fig. 5.1.2-2,3) consists 
of a fixed detector and a movable one. The movable detector 
can be rotated through an angle of 180 degrees from 90 to 
270 degrees relative to the fixed detector. Danger of da- 
maging the lead shield of the scintillation crystals pre- 
cludes more acute angles. Provision is also made for vary- 
ing the detector solid angle subtended at the source in thir- 
teen integral steps. This varies the geometrical half angle 
of the subtended cone from about 19% degrees to about 5 de~ 
grees. Due to scattering effects and the finite size of the 
source the effective half angles are somewhat broader than 
the figures quoted (See Section 5.3.3). 

The detectors themselves are on mounts which hold the 
axis 6"7/16 off the base plate. Lead shields are provided 
which interpose a minimum of 1# centimeters of lead between 
the crystals when the detectors are at right angles. The 
crystals are discussed in connection with the photomultiplier 


tubes in Section 5.1.35. 


Section 5.1.5 — Apparatus, Electronic 
General 
The electronics have four primary tasks, viz.e, detec 


tion of the radiation, energy selection in each channel, 
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time comparison between channels, and counting. 

Detection is performed by the crystal-photomultiplier 
combination and associated components and amplifiers. Ener- 
gy selection is performed by the differential discriminators, 
time comparison is performed in the fast double coincidence 
circuit. The final selection of those pulses which have pas- 
sed both the "time" and “energy” tests is performed in the 
slow triple coincidence circuit. Counting and presentation 
are done by the scalers. 

It is found necessary to separate the functions of time and 
energy selection. The reason for this is illustrated in Fig. 
Seled=-1l. It is based on the principle of operation of the dif- 
ferential discriminators which will be discussed in more detail 
later in this section. But for the purposes of this explanation 
it is merely necessary to know that one must use a finite chan 
nel width AE, i.e., the only pulses which will “operate” the 
differential discriminator are those from the linear amplifier 
in a range, E) to Eo (See Fig. 5.1.3-1). Now assume that photo- 
peak pulses arising from simultaneous gammas of identical energy 
may give rise to a pulse from the linear amplifier anywhere 
within AE; this is due to the statistical spread in the photo-= 
peak pulse height introduced by the photomultiplier tube. The 
time of firing of the "baseline" Schmitt circuit determines the 
time of the output pulse of the differential discriminator, 
and this circuit fires when the voltage of the input pulse rises 
to Ey}. It can be seen directly from the figure that a time un- 
certainty At is thus introduced in the output pulse of the dif- 
ferential discriminator. (We-normally operated with AE 45 v 
which gave rise to At = 4x107" sec.). One is thus constrained 


to use a coincidence circuit operated by pulses generated when 
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the pulse height reaches some lower value such as E,. It is clear 
that A‘t is much smaller thanAt. For this reason we use the pul- 
ses from a fast amplifier, which saturates on the first few volts 
of the photomultiplier pulse, to operate the fast double coinci- 
dence circuit. 

Block diagrams of the entire apparatus with the interconnec- 
tions for various types of runs and checks are shown in Section 


eel. 


Photomultipliers, Crystals, and Associated Components 

DuMont Type 6292 photomultiplier tubes were used. They were 
operated at about 127 volts per stage, adding up to a total high 
voltage of 1800 volts. Under these conditions exceptionally large 
output pulses were obtained (40.volts for the 722 KEV photopeak) 
with extremely good energy resolution (~6%). An ordinary cath- 
ode follower would not pass such a large negative pulse with a 
short rise time down even a moderate length of coaxial cable. 
This is Because one needs a relatively large cathode resistor 
to maintain the quiescent voltage at the cathode greater tham the 
pulse height to “be passed, and this large resistance together 
with the capacitance of the cable prevent the passage of short 
rise time pulses. For these reasons a "bootstrap" circuit was 
used, in which the cathode resistance is effectively replaced by 
a variable resistance in the form of a vacuum tube (Fig. 5.1.3=<2). 
The pulses had a rise time of about 4x07” seconds and a decay 
time of about 4,5x107° seconds. The decay time constant was 
set by the 68K resistor at the collector of the photomultiplier 
together with the total capacitance to ground at this point 
(about 30 uuf)). 

Relatively small scintillation crystals were used. They 


were 1 inch in diameter by # inch thick Harshaw Type X4L2: they 
2S gh = 





Photomultiplier Voltage Divider Unit and Bootstrap 
Fig. 8.12542 


Notes 1K-2 
1. All resistors s watt unless q 
otherwise specified. 
2. All resistors in HV divider 
plus-minus 5%. 


3. "m" = megohms 
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"OUTPUT" 
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were mounted in a hermetically sealed can made of .032" 2-98 alu- 
minum which was internally coated with MgO. Ordinary white vase- 
line was used to provide optical contact between the window of 
the crystal can and the photocathode. Crystals were NaI(T1), 

The 460 volt battery applied to the second dynode (Fig. 
5.1.3=2) had the purpose of maintaining the early stage electron gain 
both large and constant. This tends to increase the energy reso- 
lution as well as reduce the effects of high voltage fluctuations. 
The good energy resolution of this cry stal-photomultiplier combi- 
nation is illustrated in the observed partial spectrum of the gam- 
ma rays in the decay of Int14 (Fig. 5.1.33). 

There is apparently a transit time uncertainty of about 
ex10"° seconds in these tubes under these conditions, which makes 
this set-up undesirable for lifetime measurements by delayed coine 
cidences. This point is discussed more fully later on in this 
section in the paragraph on the fast double coincidence circuit. 

Atomic Instrument Company high voltage power supplies (Model 
316) were used and performed satisfactorily. 

Pulse Amplifiers. 

The pulse amplifiers (Fig. 5.1.3=4) were designed to provide 
a moderately fast pulse output together with a linear pulse oute 
put in the same chassis. In the case of the linear output the 
term “amplifier” is rather a misnomer: the term "inverter" would 
be more appropriate. The already large negative pulses from the 
photomultiplier are inverted and amplified slightly and then fed 
to the differential discriminator which requires positive pulses. 
The linear amplifier consists of a heavily fed back single 6AG7 
stage whose grid operates about 50 volts above ground. Gain from 
-1.6 to ~80 is adjusted by varying the amount of unbypassed cathode 


resistance. The pulses from the sodium iodide = photomultiplier 
2 6 es 
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combination (~40 volts, 4x10~° seconds rise time)’ when amplified 
in the fast amplifier have a rise time of about 5x1078 seconds and 


a height of about 8 volts when clipped with 4 meters of RG-7/U. 


Differential Discriminators 

These operate on the well known circuit of two integral dis- 
criminators feeding into an anti-coincidence circuit. Modifica- 
tions were made to a standard type of differential discriminator 
which resulted in the saving of six tubes and in ruggedness and 
reliability of operation. These modifications center around the 
anti-coincidence circuit (Fig. 5.1.3-5).. Briefly, the operation 
is as follows. The square wave generated by the “Base Line” 
Schmitt circuit is differentiated strongly and is fed through a 
delay line into one input of the anti-coincidence circuit. The 
delay is necessary to give the "Delta E" Schmitt circuit a chance 
to fire IF IT I8 GOING TO FIRE. Here we have two possibilities: 
(1) the "Delta E" Schmitt IS going to fire, i.e., the input pulse 
is too big, or (2) the "Delta E" Schmitt IS NOT going to fire, 
i.e., the input pulse is within the “window”. In the first case 
the square wave generated by the "Delta E* Schmitt is fed into 
the anti-coincidence circuit and prevents the delayed, differentia- 
ted pulse from the “Base Line" Schmitt from passing through, i.e., 
the "Delta E” Schmitt gates the anti-coincidence circuit off. 
Thus, in case (1) we get no output pulse. In the second case there 
is no gating pulse from the "Delta E"® Schmitt and our delayed, dif- 
ferentiated pulse from the "Base Line” Schmitt is passed through 
the anti=coincidence circuit, then amplified, and then fed as a 
trigger pulse to a blocking oscillator whose output is the output 
of the differential discriminator. 
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Fast Double Coincidence Circuit 

The heart of this piece of apparatus is the non-linear ele- 
ment, which in this case is a 6BN6 vacuum tube. The virtues of 
this particular tube are that the plate characteristics of two 
of its grids, the quadrature and limiter grids, are practically 
identical. A negative voltage of 34 volts on either grid with 
a plate voltage of 60 volts is sufficient to prevent plate current. 
The tube is operated normally with both grids in the cut-off con- 
dition ( about -5 volts ). Semi-standard positive pulses of about 
8 volts from the fast pulse amplifier in each channel are fed in, 
one to the quadrature and one to the limiter grid (Fig. 5.1.3-6). 
If they are “simultaneous” the tube conducts and the pulse is am- 
plified in a 6GAH6 stage and fed out through a 6J6 cathode follower 
stage. The fast double coincidence circuit determines the overall 
resolving time of the apparatus as a whole. The resolving time 
of this circuit can be adjusted in two ways as follows: (1) by 
“shortening” the input pulses, i.e., by using a faster, more 
strongly clipped input pulse, and (2) by adjusting the bias level 
on the grids of the 6BN6. Both methods were used: method (2) 
turned out to be extremely sensitive to the grid bias adjustment 
and made the use of a precision potentiometer necessary in the 
cathode circuit of the 6BN6. Resolving time curves obtained by 
delayed coincidences using the single channel method (Fig. 5.2.1-3)) 
and the double channel method (Fig. 5.2.1-4)’ are shown in Figs. 
Sel.d-7,8. Since the annihilation radiation has zero lifetime, 
the broadening of the double detector curve is introduced by the 
apparatus. That this is not a “channel width" effect is shown 
by comparing Figs. 5.1.3=—7 and 5.1.3=<8 in which two widely dif-= 


ferent channel widths were used and which resulted in approximately 
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the same broadening. We thus look to the photomultiplier-crystal 
combination as the source of the "jitter" which causes the broad- 
ening. That this jitter is not caused by the uncertainty in the 
time of creatiom of the first photoelectron at the photocathode 
can be shown by applying formulae due to Bell, et. al. (1B52). 

That this time uncertainty is not caused by varying trajectory 
lengths for the avalanche electrons has been shown by measure- 
ments by Alford (unpublished) using a pulsed light source. However, 
as Alford has shown again, the average transit time is a very de- 
finite function of which part of the photocathode is illuminated; 


“5 seconds longer for photo- 


the average transit time is about 10 
electrons emitted from the circumference of the photocathode than 
for those emitted from the center in the DuMont 6292's. If our 
photocathode were illuminated by several intense bursts of light 
randomly spaced in area over the photocathode, then our effect 


of broadening might be explained on this basis, 


Slow Triple Coincidence Circuit 

This circuit acts as a sort of recombining point in the flow 
of pulses through the whole apparatus, where the two functions, 
viz, energy selection and time comparison, which were previously 
separated, are rejoined (Fig. 5.2.1-1). Each of the three input 
pulses, one from each of the differential discriminators and one 
from the fast double coincidence circuit, triggers a square wave 
generator whose output pulse length is adjustable from 0.2 to 
1.5 micro-seconds (Fig. 5.1.3-9). The pulses from each of these 
three circuits are then fed to a non-linear element consisting 
of crystal diodes (1N56) in parallel. If the three pulses over- 
lap in time, the common plates of the crystal diodes rise sharply, 
and this pulse is then fed through a class © amplifier. The oute 
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put of this circuit, which is “energy selected coincidences", 
then goes to a scaler where the pulses are counted and recorded. 
The square wave generator is a plate triggered univibrator;: the 
cathode follower coupling of the two halves of the univibrator 
results in a rise time of about 5x1078 seconds for these square 


pulses. 


o 40 o 





Part 5.2 = Experimental Procedures 
Section 5.2.1 - Runs and Checks 

The various types of runs and checks to gather information 
of the source or the equipment are best understood by referring 


to the inter-connection diagrams which follow. 


Type Runs Desig- Channel Var- Fig. 
nation Width lable 5.2.le 

Total Coincidence Counting TO Wide 0 1 

Rate 

Chance Coincidence Counting Co Wide S) 2 

Rate 

Resolving Time, Single Chan—- RT=-A Wide Delay 3 

nel Method 

Resolving Time, Double Chan- RT-B Wide Delay 4 

nel Method 


The checks include operational and characteristic checks. 
Operational checks, or “alignment checks", merely serve to tell 
if the apparatus is properly tuned. These were performed daily 
before each run and occaisonally at the end of a run when there 
was reason to believe that something had "sone wrong", Characte= 
ristic checks are tests to ascertain that we are actually mea- 


suring what we think we are measuring. 


Operational Checks 

Window Set - This involved setting the photopeak pulse heights 
of the 722 KEV or the 556 KEV gamma rays in each channel to a pre- 
determined height, then by means of the pulser to set the lower 
window edge at a pulse height corresponding to 430 KEV. This fi- 
gure was arrived at from the Compton scattering chart of the gam—- 


ma rays present in the decay of tnt}* 


, Fig. 5.2.1-83 see also $ec- 
tion 5.5.5. Then, using the pulser again, delta-E was adjusted 
to a predetermined value to include both lines (722-556 KEV) in 


h ch e 
each channel .& : 





Explanation of Abbreviations Used In Figs. 5.2.1-1,2,3,4,5,6,7 


PM 
CF 


PMHV 
(amC’) 


LO 
FO 


DIFF. 
DISC. 


FDCC 
STCC 


SCOPE 


Photomultiplier Tube 
Cathode Follower ("Bootstrap") 


Photomultiplier High Voltage Power Supply 
(Atomic Instrument Company, Model 316) 


Linear Output of Pulse Amplifier 
Fast Output of Pulse Amplifier 


Differential Discriminator 


Fast Double Coincidence Circuit 
Slow Triple Coincidence Circuit 


Scaler 


Cathode Ray Oscilloscope, Tektronix Model S1lA 


All interconnection lines represent coaxial cable RG-7/U unless 


otherwise stated. 
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Slow Triple Coincidence Circuit Cout-Out Test - This tests 
the overall performance of the energy and time discrimination func- 
tions of the apparatus. The input to each differential discrimi- 
nator and each input to the fast double coincidence circuit are 
removed in turn. The output of the slow triple coincidence cir- 


cuit should disappear with each removal. 


Scaler Counting Efficiency Test (Fig. 5.2.1-5) - With the 
pulser set at the normal single channel frequency and so that 
the pulse height falls within each window, the scalers should 
count together. Two scalers were compared at a time, viz., 
$2:83, $2:84, $3:84. It was found that they counted together 
accurately to within plus or minus .01%. This test also par- 
tially checks the performance of the differential discriminators. 


Characteristic Checks 

Cable Loss (Figs. 5.2.1-6,7) -— Attenuation in the 50 meter 
RG@7/U delay cable could cause the measured chance coincidence 
rate to be less than it actually is. Using the set-up shown in 
Pigs O.2.1-6, there is an observed dimunition of counting rate 
as the delay is increased, whereas no such droop is observed in 
the energy sensitive case as shown in Fig. 5.2.1-7. The former 
effect is evidently due to attenuation of the lower energy pulses 
in the fast amplifier output which are normally just above threshold 
for the fast double coincidence circuit, and which do not trigger 
the differential discriminators in any case; whereas those pulses 
which are energetic enough to fall within the window are still 
above the fast double coincidence circuit threshold in the latter 
effect. We conclude that cable loss in the energy sensitive case 
is negligible. 

sn OES 





Scattering - Two ways in which the scattering of gamma rays 
could adversely affect us are (1) by scattering from crystal to 
crystal which would cause a spurious coincidence to be registered 
under certain conditions, and (2) by a change in the scattering 
conditions in or near the source which would then be unaccounted 
for by the geometrical corrections. (1) is accounted for in two 
ways: (a) by the detector lead shields (Section 5.1.2), and (b) 
by the location of the window base line. (a) is useful at 90° 
but of no use at all at 180°, whereas (b) is most effective at 


180°, Fig. 5.2.1-8 shows the enrgy of the scattered quantum versus 


scattering angle for the gamma rays present in the decay of tnil4, 
For a crystal to crystal scattered coincidence, the gamma ray must 
scatter at least through 90°, The window base line of 430 KEV is 
chosen so that in the worst possible case (the weak 1300 KEV com- 
ponent), the window excludes all rays scattered through more than 
78°. 

The effective broadening of the solid angle of the detectors 
due to scattering from the shields, in the source, etc., is ano-~ 
ther matter and is discussed under "Geometrical Corrections", Sec- 
tion 5.3.3. However, these corrections would not correct for (2) 
above, i.e., VARIATIONS in the scattering conditions due to frost 
on the source holder, variation in the phase of the source, and 
the presence of a thermo-couple in the source, etc. In order to 
check for these variations in scattering conditions, two angu- 
lar resolution runs were made using Na&<, one under conditions of 
maximum scattering, i.e., using a thermo-couple in a heavily fros- 
ted source holder, and one under minimum scattering conditions, 
i1.e., at room temperature with no thermo-couple. Effective increase 
in scattering conditions would show up as a broadening of the 


angular resolution curve. Results, however, showed no detectibdles 
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difference in the curves. We conclude, therefore, that the ef- 
fective scattering conditions remained constant throughout the 


experiment. 


Section 5.2.2 -—- Sources 


114 and Na22. Both of these iso- 


Sources were prepared of In 
topes came from Oak Ridge in the form of chlorides in dilute aqueous 
solutions of hydrochloric acid. 

ee came as Intl. ima .38N solution of HCl. Two separate 
shipments of 10 mc each were received, (1) having a specific acti- 
vity of 171 mc/gm and a concentration of 1.90 mc/ml and (2) having 
a concentration of 87 mc/gm and a concentration of .96 mc/ml. 


114 sources were prepared, three of them of 


Four separate In 
the extended source geometry shown in Fig. 5.2.2-1 and one (200 
uc) merely evaporated to dryness then sandwiched between to "1/8 
aluminum plates for lifetime measurements using another apparatus. 
Of the three sources using the “extended geometry”, two were in 
liquid form and one was in solid form. The solid was obtained 
by evaporating 1 me of shipment (1) to dryness in the "solid" type 
extended geometry. lucite vial (Fig. 5.2.2-1). Of the liquid 
sources, one was of "1 centipoise” viscosity, i.e., used just 
the way it arrived from Oak Ridge (shipment (1)), and the other 
was of temperature-variable viscosity, prepared from shipment 
(2) as follows. Two ce of the stock solution were evaporated to 
~O5cc, then 6P pure glycerin (greater than 99% pure) was added to 
make 1.95 ec. The remaining volume of .05 ce was filled by the 
thermo-couple. This solution was duplicated using non-radioactive 
indium for experimental determination of viscosity versus tempera- 
ture. Relative viscosities were measured using a simple "rate- 
of-fall" viscosimeter. The results are shown in Fig. 5.2.2=2, 
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Points for zero degrees and above were from (1M53); points below 
zero were experimentally determined. Fig. 5.2.2-3 shows a photo- 


114 _ Glycerin 


graph of the source mounting arrangement for the In 
source. The activity of each of the liquid Int14 sources was 
about 1.5 mc. 


“2 sources were prepared, each of about 


Three separate Na 
100 uc activity. Two of these were prepared in “extended source" 
containers, and one in a “point source" container (Fig. 5.2.2-1). 
Of the two in "extended source” containers, one had a thermo-couple 
and one did not. This was to check on variations in scattering 
conditions as discussed in Section 5.2.1. Fig. 5.2.2-4 shows a 


ee point 


photograph of the source mounting arrangement for the Na 
sources 

It is worthwhile at this time to mention something about the 
schedule followed relative to the sources. First of all, the com- 
plete, undisturbed correlation was measured, using the dilute 
aqueous solution. Following this, the various disturbed "aniso- 
tropy only" measurements were taken, interleaving from day to 
day ina random manner the source to be measured. Also, the "ani- 
sotropy only" measurements for the undisturbed correlation was 
measured about midway through the disturbed "anisotropy only” 
runse In this way it was hoped to minimize the possibility that 


long term "drifts" in the electronics could be mistaken as an 


effect of the source. 
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Part 5.5 - Data Reduction 
Section 5.3.1 - Errors in General 

The raw data as it is measured has in it many errors. We 
can divide these errors into five groups: 

1. Systematic Errors - Errors due to malfunctioning or misa- 
lignment of the apparatus, such as electronic drifts or source 
eccentricities. 

2. Inherent Errors ~- Errors inherent in the design of the 
apparatus, present even when adjustments are perfect and opera- 
tion stable, such as errors due to finite detector and source 
size. 

So. Statistical Errors - Errors observed due to the random 
nature of the events counted. 

4, Lifetime Errors - Errors due to the continually changing 
strength of the source due to its radioactive decay. 

5. Decay Scheme Errors - Errors in the observed results 
due to the nature of the decay scheme, such as the presence of 
an annihilation coincidence rate at 180° due to the presence of 


@ positron component in the decay. 


Corrections must be made for each of the above errors. The 


sections in which they are discussed are listed below: 


Error No. Section No. Title 

1, 4 Dedek Data Normalization 

“a Deeded Geometrical Corrections 
O Debe4 Statistics 

4 Saaer Decay Scheme Corrections 
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Section 5.5.2 -— Data Normalization 

Assume that one has a source which decays in a single gamma- 
gamma cascade. Let the source strength be No decays per second. 
One observes a coincidence rate as a function of the angle between 
the detectors. Assuming no coincidences arising from crystal- 
to-crystal scattering (Section 5.2.1), the observed coincidence 
rate is the sum of two independent rates:: 

(1) a true coincidence rate, N,., due to the two gammas in 
the same decay, and 

(2) a chance coincidence rate, N.,, arising completely by 
chance from separate but simultaneous decays. 

Letting the observed, or total, coincidence rate be N_.(0), 
we have: 
oy se (ea) ="8,(0) 4°. 

Furthermore, considering also the time dependence, 
Py wh (ot) Ny (4 0G )D keheuet Nee Nahey 

where W(C) is some function of the angle between the de- 
tectors, \ij is the probability that an isotropically emitted 
gamma ray, 4, will be detected and recorded in channelJ, and 
XY is a further reduction factor due to the electronics. It in- 
cludes the efficiencies of the fast double and slow triple coinci- 
dence circuits, and the scaler counting efficiency. 

The chance rate is expressed by, 
(3) Noo(t) = 22n5(t)C dat deall Mbt Azb] X 

where & is the “resolving time" of the apparatus, i.e., the 
maximum amount of time by which two pulses can be separated and 
still considered “in coincidence”. 

If we express the single channel rate by N¢i , then, 


(4) Ni) = Not) ee tA2E I i= a,b 


ay ie 





Combining (1), (2), (3), and (4) we come down to: 


(5) N to (0; t) — W (6) hice ly aT Cy 
Na(t)N p(t) No (t) Oye Peale ye ual 


But we are interested in the true coincidence rate. By a 


slight transformation (5) becomes: 


(6) Nre (et) = Nao (Ost) — New (t) = We) Nott) Nott) | (dadaist haadae) ¥ 
Neo (t) (He: \eah | 


Now, to a first approximation the factor on the right in 
brackets is independent of systematic errors consisting of small 
eccentricities of the source and small electronic drifts. This 
will be shown later in more detail. 

The ratio of W(O ) at two different angles O, and G2, at two 
different times bl. and bocan be extracted from (6) by dividing the 
true coincidence rates at these angles, remembering that the fac- 


tor in brackets is to first order a constant. We then have: 


(7) W (62) — Nre (B., ity) Nal Os, b:) Nol(@.ty No (t2) 
W (A,) Nre (O81, ti) Na (62)t2) N4(2,t2) No (68) 


which is our desired form, normalized against the small errors 
already mentioned. It is also the form in which we are interested 
in the correlation function, W(0), ices, as a ratio rather than 
as an absolute value. 

In order to show that the factor in brackets in (6) is approxi- 
mately constant we begin by defining, 


Wee 
a,b 





(8) (i, = aa i,k 
pli, , es ; 


and examining Nia little more closely: 
(9) Nit = Usai E4i Kit 

where (MO is the effective relative solid angle presented by 
the detector to the gamma ray, and 


€ is the detection efficiency of the detector, or counter 


aw AO na 





as it is often called, for the gamma ray, and 
o is the channel factor, i.e., the fraction of the detected 
gamma rays which fall within the window. 


Then, 
a ie ge 


(10) Misi _ GIse €s0eoEL are: 


Wki Exc Ske 

Now, the assertion is that (10) is insensitive to small elec- 
tronic drifts and small eccentricities in the source. For although 
both Sit and WJi will change with the fluctuation in various vol- 
tage levels, the ratio ae and at will not change if the win- 
dow is properly chosen. Similarly, wiéiis a function of O for an 
eccentric source, but wi (0) is a constant for small eccentrici-~- 

ki CO) 

ties. 

By making use of (8) we can express the factor in brackets 
in (6) as follows: 


| ( \Losz1 + pb, 21) ¥ 
(1+ pasa \(1+ pib,ou ) 


From the discussion in the previous two paragraphs, we see 
that the {vs are approximately time and angle independent. Assu- 
ming an approximately constant coincidence circuit efficiency, 
expressed by ¥ , we have thus shown that the factor above is 


approximately time and angle independents; hence (7) is justified. 


= Se 





Section 5.5.35 —- Geometrical Corrections 
General 

The directional correlation function, Wig , as derived by 
Hamilton (1H40), represents a microscopic phenomenon, viz., the 
relative directional probability for the decay of a single nucleus, 
The apparatus, however, observes a “macroscopic phenomenon macro- 
scopically", i.e., it observes a finite volume of many millions 
of nucleii, whose exact location is uncertain to the extent of 
the linear dimensions of the source volume, and it observes these 
with extended detectors which introduce an additional angular un-= 
certainty. Hence we are not precisely sure where the gamma ray 
we count came from nor exactly where it went. What we measure 
is a weighted average of W(q\ integrated over the detectors and 
through the volume of the source. We wish to determine how the 
coefficients of our measured function, W@),are related to the 
coefficients of the “microscopic” function, W(9), where 4 is 
the angle between the axes of the detectors and, as usual O is 
the angle between the propagation vectors of the two radiations 
involved in the cascade, These relations we call "geometrical 
corrections", ~~ 

The reader is referred to Appendix III for a general treate- 
ment of the philosophy of these corrections. There we present 
the basic arguments for our consideration of the geometrical cor= 
rections under the following two headings:: 

(a) Finite Detector - Point Source, and 
(b) Point Detector = Finite Source. 

In this section we will discuss (a) in some detail, but 

merely quote the results pertaining to (b) which is treated more 


extensively in Appendix III. 
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Finite Detector - Point Source 

In 1F51 Frankel shows that the effect of the finite detector 
solid angle subtended at a point source is to cause a simple reduce 
tion in the magnitude of the correlation coefficients. That is, 
the observed function, 

(1) W/o) = px OQ ox HM oq Pz, (cose) 

where the A, are the coefficients of the “microscopic" 
function, 

(2) Wig) = oan A2% P22, (cose) 

In 1652 Church and Kraushaar indicate a convenient method 
of experimentally determining the Q 5, using the “combined” re- 
lative angular resolution curve, H (6) » which is the directional 
correlation of annihilation radiation. That is, 

(3)) Q5. = \ H (9) P +x (coseid (case) 

where H(g\ is the “combined" relative angular resolution 
curve using annihilation radiation, and © is the angle between 
the axes of the detectors. 

In 1L53 Lawson and Frauenfelder point out essentially that 
the H (6) of 1052* does not give a true picture of the overall re-= 
sponse under most conditions. This is because the individual 
detector resolution curves, which Lawson and Frauenfelder call 
Evi} , X being the angle between the detector axis and an 
element of detector solid angle, dw », are strongly dependent on 
gamma ray energy, Eg , and the window base line energy, Eh, 
i.e., the minimum gamma ray energy accepted by the pulse height 
selector. In other words, gamma rays of different energies "see" 
apparently different solid angles, and in addition any change in 
“base line” discriminator setting will also cause apparent solid 
angles to change, All this points up serious limitations to the 
annihilation radiation methods of I€52. However, calibrating 


. £2 < 





each detector with gamma rays of energies identical with, or close 
to, the ones in question, as required by the method of 1L53, is 

by no means easy;' it is sometimes impossible to obtain a single 
gamma ray of the required energy and effective collimation is dif- 
iecult. 

The purpose of this discussion is to point out a refinement 
of the method of Church and Kraushaar which will remove some of 
the difficulties indicated by Lawson and Frauenfelder, and still 
retain the experimental simplicity of the former method. 

The proposed method is based on the assumption that we can 
imagine a combined angular resolution curve, H’ , based on hypo- 
thetical gamma rays of the energies in the cascade but which 
have an agular correlation like the 511 KEV annihilation radia- 
tion, i.e., a delta function, and that aside from an arbitrary 
normalizing factor, 

(4) fio 4 <., cama fe Hoo, Eg3,© tn, E bz) 

where Eq, and Eqa2 are the energies of the cascade gamma 
rays and E43 is the energy of annihilation radiation, 511 KEV. 
Eh, is the window base line energy used in practice and Eb 
and B4, are "calibrating" window base line energies chosen so 
that the combination: 

Bat? EB wt "sees" the same solid angle as Bast? E iy ; 
and Eo,+> Ey "sees" the same solid angle as Eqs E42 ° 


Or, expressed symbollically: 
(5) (Cea (2 aes Os (Eas, Ebr) 
and 
Ba (Een, Ee) ~ O3(Eqs, Eb) 


where Bc is the angle of Compton scattering which degrades 
the gamma ray from Equ to Ey, © C=1,2 These calibrating 
base line energies, Bb, and AY are determined using Fig. 


Yecel=8, For reasons discussed in Section 5.2.1, a value of Eh| 
oe 1a 





of 450 KEV was used throughout the experiment, 

In other words (4) and (5) assume that the half width of the 
resolution curve, but not its essential shape, is different for 
various Eo,» Eh combinations. This assumption seems justified 
up to about Eq =l MEV. There is also the implicit assumption 
that the differential cross section per unit solid angle for Comp-~ 
ton scattering remains approximately the same; 

(6) oF (ai) pe £0 (Bs) Eales 

Thus, the procedure is to determine from Fig. 5.2.1=8 
the base line energies Ou and BAe which bear to E 93 the same 
relation (via the Compton scattering mechanism) as does E4, to 
Eg, and Ea, (5), then measure HQ,Eq3, Et, E02) and by (4) use it 
for H(a,€91,Eq., Eb.) » Obviously the method is only an approximation, 
but the smaller the difference E,,-Eqx , the better the approxi- 
mation. C =1/1,2 

Now let us apply (3) with a slight modification as follows: 

Since we measure coincidences in two ways indicated in the 


following table, both of the following combinations apply: 


Combination No. Gamma Ray No. Detector 
_ iL a 
z db 
z 1 db 
ZL Q 


We must take two “combined” angular resolution curves, 
H-(0,£q3,E'5,Eb2) » Using annihilation radiation where < represents 
the combination number above, and the respective window base lines 


are given belows 


—_—al - a ; = _ 
—- a= = a a 
ew ~ 
nm i/i= = i. 
a 


o 
"4% 





Combination No. Channel 1 Channel 2 
(Detector a) (Detector b) 


Base Line Base Line 
1 Bb E bz 
2 E42 BI 


Then (3) becomes: 


(7) (Q 2« =tLdiact darn] 


where: 


(8) Meek = Pek (cosa) Hi (0) d (case) LalyR 


It turns out that in this apparatus that: 
(9) Hi(o) = H2(6) 
Hence by (7): 


(10)! Q 2k = i Pak (c.as@) Hi (@)d (cas @) 


We thus have a method for calculating the coefficient correc- 
tion factors of our "Finite Detector = Point Source" correction. 
The method retains the experimental simplicity of 1652 while ree 
moving to first order at least the major objection raised to 10652 
by 1L53. 


Point Detector - Finite Source 

As remarked previously, we shall only quote here the results 
which are considered in more detail in Appendix III. 

We observe a function, W(@) » where ) 4s the angle between 


the detectors, and W(@) is given by: 


(any) W (y= ee Care (cos @) 
K 


(Note that « may possibly include odd integers). 
We wish to determine the relationship between the ies and 


s 
the Qa Axx Sof (1) which we will designate by Bak ° 
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Then from Appendix III: 


(12) Co = B+ § (Bo + = Bet _ Be) + eg B,- a B,- = B,) 
¢, = f ( Be- SB,- BL) + (Bot Be) 
C,= Bat § (2B,) + YC $B,-2 8B, ) 
= g, (8,- 3 8,) + ¥(4B,) 
oe) = iat §.¢3 B, ) + Yo-Z By) 
Cl = §,( SBy) 


Where & is the ratio of the radius of the source volume con- 
Sisting of a centered right circular cylinder coaxial with the gz 
axis of the coordinate system to the distance from the origin to 
the detectors and J, is the ratio of the “half-height" of the cylin- 
der to the origin-to-detector distance. Actually, the Br's appea- 
ring in (12) do not yet correspond exactly to the QauAzsof (1), 
but when normalized to Ro={ they will, It is almost too obvious 
to mention that the finite source does cause “mixing” of the coef= 
ficients, and it thus behooves one to keep § and §, as small 


as the requirements of the experiment will allow. 


Other Geotte trical Sources of Error 

& pitfall that must be guarded against is the distorting 
effect of the lips of the lead shields of the detectors, which 
have the effect of making the crystal response more sensitive 
to the exact position of the source than would be the case in 
the absence of the lead. (For this reason there would have 
been some advantage in the specific case of eqit* where the gam— 
mas being investigated are fairly close in energy, in eliminating 
the lead shields and in relying upon the available energy dis- 


crimination to avoid crystal-to-crystal scattered coincidences 


(Section 5.2.1), but in general this method fails when the gammas 
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are far apart in energy). 
The distorting effects of the lead shields take place in two 
ways: 
(a) The lips of the shields accentuate the spurious corre= 
lations arising from source eccentri- 


cities in the plane of the detector 





Lean woe 
axes, as indicated in the figure at SHIELDS 
right. For an unshielded detector 
the principal effect of this horizon- 


tal eccentricity is that caused by 


variation in source-detector dis-— 





tance as the angle between the de~ 
tector axes is varied; with lead 
shields this effect is accentuated by the undesired shielding of 
the detector by the lips of the shield at certain angles. 

(bd) For a point source which is off center in the vertical 
direction, the lips of the shields, by shielding the detector in 
an asymmetric manner relative to the source, cause the effective 
detector angular resolution curve to lose its symmetry with res- 
pect to the detector axis. The fact that there is no mixing of 
coefficients in the "Finite Detector = Point Source" correction 
(1), is based on the detector axial symmetry of this resolution 
curves hence in this case if we apply (1), we will introduce er- 
ror, in this case not only does the effective resolution curve 
lose its aforementioned axial symmetry, but its half width changes 
also as the source position is varied in the region where the op- 
tical view of the detector from the source is masked by the lips 
of the shields. We have spoken so far of a point source, but what 
we have said applies equally well to each element of a finite 
SOUrCE. 
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Of the two types of error, (b) is the more insidious in that 
it has no well defined symptom which arises in practice. Type (a), 
for instance, manifests itself immediately by an anisotropy in the 
single channel rates, whereas (b) shows up only as an error in the 
results. 

The precautions taken against these errors are careful mecha=- 
nical and electronic adjustments. In case (a), we center our source 
in the plane of the detector axes by observing the single channel 
rates as a function of detector angle and choosing that position 
which yields the most isotropic single channel response. In the 


ee point source was aligned 


vertical direction, case (b), the Na 
at the position of maximum coincidence counting rate, which is 
very sensitive to vertical eccentricities, at a relative detector 
angle of 180 degrees. The same procedure was followed for the 
extended Na source, and the position thus found was used for 
the Inité4 sources which used the same geometry as the extended 
ie” SOUPrCe. 

One might think that errors arising from case (b) might best 
be handled by measuring the resolution curve associated with the 
a extended source and using this in a calculation as per (1) - 
(10) above, but from (12) we see that the coefficient mixings 
introduced by the source extent would cause error if we used in 
(10) the results obtained with an extended source. The proce- 
dure followed was first to take a combined angular resolution 


rife 


curve, as defined by (3) and (9) above, using the Na~~ point 


source, then determine the "Finite Detector = Peint Source" cor= 
rections from (10) above. Then the coefficients of the tnil4 
correlation, Bo, , uncorrected for finite detector size, were 


extracted from the observed coefficients, C, , using (12). The 
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true coefficients, A... , were then determined from the B2, by the 
application of (10). In applying (12), there is an ambiguity in 

o> corresponding to the fact that the "half—height" of the source 
is defined in part by the aperture in the lead shield of the de- 
tector: but this correction was small enough not to be sensitive 

to assumptions on this point. The same calculations also assume 
that the "Finite Detector = Point Source” corrections, based on 

the detector angular resolution curves, were the same for source 
elements at all heights. That this is not exactly true is discus- 
sed under case (b) above, but since we have an extended source, 

for small displacements of the source center, there is always a 
source region of the same size which has a "view" of the detectors 
uninterrupted by the lips of the shields. Thus, the resolution 
curves from source elements in this region, which contribute much 
more to the coincidence rate than those elements in the "shadow" 
zone, still retain their symmetry and shape necessary to the ap-= 
plication of (10) as previously mentioned. We find then that with 
careful vertical positioning of an extended source, the corrections 
are stationary with respect to vertical motion of the source center 
and thus the errors arising from case (b) are second order in na- 


ture. 
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Section 5.3.4 = Statistics 
By the nature of the processes involved the measurements we 
take are samples from a probability distribution. The correlation 
function itself, Wea), is a probability function, i.e., it repre- 
sents the relative probability for the emission of two successive 
quanta at an included angle, 0 , between their propagation vece 
tors. if, for example, we take two successive readings of the 
anisotropy, A, where A = W(180)/W(90) = 1, there is but a vanishe 
ingly small chance that they would read exactly the same, but this 
does not mean something is wrong with the equipment. It means 
that we are sampling from a probability distribution of A, cene 
tered we think at some mean value, and if we take enough samples 
we think we can specify that mean with a given degree of accuracy. 
In connection with statistical accuracy, there are several questions 
that may be asked. Among these are?’ 
(1)’ Is there an optimum source strength? 
(2) How long must one take data to obtain a given statisti- 
cal accuracy? 
(3) Is there an optimum division of time to be spent on the 
measurement of the chance rate? 
Notation 
Let Nc =Number of counts in time, Tse. 
Let Mc = Number of counts in unit time. 
Where ifs: 
c= TO We mean TOTAL counts. 
(=> TR We mean TRUE coincidence counts. 
L= CH We mean CHANCE coincidence counts. 
(= © We mean SOURCE STRENGTH in Number of Decays. 
First of all: 
(1) ies sine Ste 
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The apparatus measures the number of TOTAL coincidences. 
This number is made up of two components, (1) the TRUE coincidences 
due to cascade radiations proceeding from the same nucleus, and 
(2) the CHANCE coincidences due to radiations arising in separate 
but simultaneously decaying nucleii. Then we can write: 
(2) Nro = Nra + New 

and 
(3) Wro = Wrer li Mow 

If we define (ANre) a8 the variance, or mean square deviation 


in Nra, then it is desired to see the functional dependence of 





the fraction a » which we define as d 
T 
Since the variance in the TOTAL and CHANCE coincidences are 


independent:: 


(4) DONG as y (A N-ra) JL ANagle 


If we assume a normal distribution for our directly measured 
quantities, then we have from elementary statistical methods:: 
(5) ANé¢ =-[ Ne C= ronrctH 

Let us make one more assumption: that we have available a 
total time, | , “which is to be divided between [+0 measuring the 


total rate and Tee measuring the chance rate. That is; 


(6) . a T +0 + Ales 


Note that: 
Pro | re 
Define: nee 





— - 
That is, & is the fraction of the total time spent on the 
measurement of the chance rate. 
Then: 
Me lr. = Ci-x) J 
And: 


(0) Ton = XT ee 





Now during that fraction of the total available time, | , 
to be spent on the measurement of the TOTAL rate, we measure the 
following number of TOTAL counts, by (1) and (9): 
(11) Nro = M0 (1-«) T 

And in accordance with (5), since this is a directly measred 
quantity it is equal to its variance, or: 
(12) (ANto)* = Neo 

We recall again from (5) that the mean deviation in the 
MEASURED chance counts would be just the square root of those 
measured during time oT ; orymenX | » but we must remember that 
we do not necessarily spend the same length of time measuring CHANCE 
as we do measuring TOTAL. Hence this mean deviation must be nor- 


malized to ({-%)T , the time spent on measuring TOTAL. Then we get: 


(13) MNca = CA) I menXT 


Or: 


De 
(24) (ANew* = EO men 


Combining (2), (4), (12), and (14) we have: 


(15) ANrr =) = + IL 
Nra ie Mra T &(1-«) 
Where: 
(16) mn = ics 
+e 


We are now in a positon to answer question (3) above. There 
is indeed an optimum division of time, for d has a minimum with 


respect to & . It turns out that the optimum & is given by: 
(17) a = nl fen -| | 


Substituting (17) into (15), we finally obtain: 
me) d = | (fitter +i) 
| fita-' —1 ]Mre T 


ay ee 





(18)) represents the obtainable statistical accuracy if time 
is used to the best advantage. On the other hand, it has no mini- 
mum with respect to-1- . However the normal way of changing Ju 
would be to change the source strenth, in the process of which 7 
and@mye are proportional to each other. A more useful form of 
(18))} in which the dependence of Mre andm onM. appears expli- 
citly will now be developed. If we now make the substitutions 
(See Section 5.3.2):: 


(19) Mey — Se "G NAS 
(20) Mre = KriW(6) Mo 
Where the K’s are constants, then by (16)): 
= Ki 
—™) S Se ere 


Futting (20) and (21) into (18) we obtain: 
(29) a — ide) oe ) 4. Ke W(6) +1 | 





Ki W(6) Ki T Me 


which has no minimum with respect to Ma except as M.->©. 
Hence we have answered question (1) above in the negative;: there 
is no optimum source strength. There is however, an effective 
upper limit on Mo - This is determined either by the counting 
rates at which the electronics cease to function efficiently or 
by safety considerations, whichever yields the lower value (See 
Section 5.2.2). 


We answer question (2) by transposing (18): 


meee se tf eet OI 
Wred* qi + 2c =| 
This function is plotted in Fig. 5.3.41 as | vs. M;, for 
various values of the parameters d and n . Thug, if we have 
a given source and resolving time, G » characterized by a given 


Mra and J. , and desire a given accuracy, d » the necessary 
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total running time may be read directly from the figure. 
We insert here a slight reformulation of (18) which is use- 
ful in that it indicates how d(n) changes as we change J. by 


whatever means? 


(24) don) ~ dco, feet pl+no +l 


fie 


The statistics we have considered so far have been those of 
a single large number. However, in this experiment our interest 
lies finally in the uncertainty of a ratio, close to unity, of 
two large numbers. What we finally want to know is the anisotropy, 
A, of the correlation where A is defined, neglecting the small cor-=- 


rections of Section 5.35.23: 


— Nradizo’ 
(25) A, = SS se 
Nre (90) | 
It can be shown that the uncertainty in A is given bys 
(26) AN = (A 4]) | A Neecrtes | 4. [Ate csen)" 
UN TRI N races 


Since the two terms under the square root are each approxi- 


ce 
mately equal to d » we have: 


(27) AA x Scale d 
-_~ 


For example, assume A = .14, what TRUE coincidence statis 
tical accuracy would be necessary to reduce the relative uncertainty 
in A to 1%? The answer from (27) is d= .000869, And using (23) 
assuming W=—1/5 and a true coincidence rate of 10 counts per 
minute, to attain this accuracy would require 528 hours of rune 
ning time! For this reason the uncertainty in A is somewhat large. 
It was the practice in this experiment to run until d= 8+ 3x10 =o 
which yields 6x10” < a < ieee dertt 


aon = 





We now wish to investigate briefly the interpretation of sta- 
tistical samples. That is, if we take some readings each of which 
has a certain plus or minus attached to it, and if we draw in the 
curve which gives the least squares fit to these readings, then 
we ask the question: “How well do my observed points ‘fit’ the 
curve drawn through them?” Or: “What is the probability that a 
sample set of readings each chosen completely by chance, having 
in mind the probability distribution at each point on the curve, 
would give no better fit?" To answer these questions we use the 
well known X™~ test, whose table answers this question direct- 
ly. If then, it turns out that, say 99 times out of a hundred 
chance alone would give a better fit, then we have a strong sus- 
picion that our apparatus is not performing its function of "ran- 
dom sampling” properly. However, it must be borne in mind that 
we are not sure that the apparatus is malfunctioning. Perhaps 
we have just measured that sample which is the one in a hundred 
instance when chance would give just the fit we measured, or worse. 
And if we take several hundred samples we are bound to run across 
this situation. This is to say that the probability distribution 
as given in the ™ ~ table is spread evenly from zero to one. 
In such a situation, if a repeated group of readings give the 
same probability, then in only one case in ten thousand would two 
such random samples appear in succession with such a Ke - then 
we would be well justified in looking to our equipment for non= 
random errors. One characteristic of the 2 a test, however, 
is that it is concerned with three or more readings, or samples. 
We develop below a way of looking at just two readings. 

Let us ask ourselves the following question: "Given a nore 


mal distribution, what is the probability, PC xo) , that two suce 


=- 65 @ 





cessive samples taken from this distribution will lie X. or 
greater standard deviations apart?" The following expression (28) 
gives the required probability, where the terms used together with 


a “picture*® of the _ ees Se is shown in the diagram. 


(28) Pex = ory Omnis | = ) Pow dx | dx 





is plotted in Fig. 5.3.4=2, and is gutionts below for 


some rational values of Xo , 


x van) Pie) 
Oo 1.000 
eLo 870 
g00 e lOO 
e lO 2595 
1,00 479 
1,50 287 
NOL, » L57 
3200 2055 


In using this data one must remember that one started with 
the assumption that the two samples were from the same normal 
distribution; it does not neccessarily follow that the inverse 
is inferable, that given two samples of a given separation, that 
Fig. 5.3.42 gives the probability that they are from the same 
distribution. It may happen that there are two separate distri- 
butions very close together, and the probability that the points 
are from the same distribution in this case would be less than 
indicated in Fig. 5.5.4-2. However, with this limitation in mind 
it seems that Fig. 5.5.42 is useful in that it gives at least 
an upper limit of the probability that two points are from the 


same distribution. 
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Section 5.3.5 - Decay Scheme Corrections 

The observed counting rate must be corrected for certain 
features of the decay scheme. In this experiment the most im- 
portant correction of this nature, and only one of any conse-=- 
quence, is that due to the annihilation radiation arising from 
the positron component in the decay of In’ <. 

The decay scheme as verified by Johns et. al. (1/54) is 
shown in Fig. 4-1. The annihilation radiation of the positron 
component although exceedingly weak is sufficient to cause a 
troublesome, unwanted coincidence counting rate at angies near 
180 degrees, which if not corrected for will yield an erroneous 
value of the anisotropy. The method by which this problem is 
attacked is discussed at length in Appendix I. The result is 
that the positron corrected anisotropy is obtained ty simply 
subtracting a factor, 6 , (6<\ ) from the observed anisotropy. 


~ 6 


Bach of the observed points are then corrected for posi- 


sti 


(1) N= A 


trons by applying the following relation: 


(2) Cog (0le= Co (8) — 8 (8) Co (401 


where’ 
(3) Ge) = H (2) 8 


where H,(@) is the positror annihilation radiation angular 
correlation function normalized so that H,(180)=1, using 
the extended source geometry and a window base line of 4350 KEV 
in each channel. 

Another possible source of error which bears mention is 
the presence of the 576 KEV transition from the 1856 KEV level 


- 67 = 





; 
in cali, However since this transition is very much weaker 


than the 722-556 KEV cascade (relative intensity = .03) (1.154), 
it is neglected. All other gamma rays presnt in the decay are 


effectively excluded by energy selection. 
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Section 5.3.6 - Curve Fitting 
After data reduction the observed function W'''(6) was fit 


to the observed points by a least squares procedure. Where: 


eee wt (Co; — adt' 4 a' cose a al! cos*6 
The parameters, ag! ; ag! 's ans are varied so as to 


minimize X~ where: 


Pips ml a 
x ml A Nte( 82) 


Since approximately equal times were spent on each point: 
Me 26h oN Nova (Qc) & ANTR(CG)) 


Then we can write: 
—- Ws — ¢ 
(4) X* ~ F | Nra(od—W (ea | 


A total of 10 points were measured, 5 each in the ares from 
90 to 180 degrees and 270 to 180 degrees. Since there was no 


phase shift in the correlation the points were numbered as follows: 


e 90,270 120,240 135,225 150,210 180,180 
i ™ a: 2 3 4 
By minimizing (4) with respect to the parameters we arrive 


at a perscription for the parameters: 


(5) ag'' = = ( Sly 4 oy meee + 34) 
alt! = <i (-108y, + 26y, + 80y, + 54y, — 5274) 
fe cogil) r =. r = ae Say, - 
_' => = ( Lem. By, 16y,, y, + l6y,) 
Where? 
(6) Ys = N+-a (82) 
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Section 5.2.7 = Data Analysis Precedures 

In this experiment we observed the following: 

(a) Complete Correlation - Enough points were observed to 
make a specific determination of the correlation coefficients. 
This was done only in the case of the undisturbed correlation 
(dilute aqueous solution of In€13). 

(b) Anisotropy Only - Measurements taken only at 90-180 
and 270-180 degree points. This was done in the case of all 
the disturbed correlations (glycerin solutions, salts, and 
frozen solutions), and also in case of undisturbed correlation. 

The procedure in case (a) is as follows: 

1. Normalize the data, i.e. correct for systematic errors 
as Giscussed in Scction 5.3.2. 

2. Obtain a least squares fit and calculate tentative 
coefficients. 

3. Using results of "2" make first approximation of po- 
sitron correction. 

4, Using results of "3" correct observed points for po- 
sitrons and recalculate the coefficients. 

5. With these new coefficients recalculate positron co- 
rection for the second approximation. 

6. Repeat “4" and "5" until two successive approximations 
of the coefficients do not change the positron correction. 

7. Gorrect for finite source size. 

8. Gorrect for finite detector angular resolution. 

9, Find correction factor for anisotrony in undisturbed 
case and apply this to results of the measurement of “anisotrory 
oniy™ in the unspoiled case. 

10. From "9" and "8" determine final coefficients by fin- 
ding the center of the areas of overlap ina plot of one cor 


relation coefficient versus the other (See Fig. 5.4.a+1). 
-_ Ve 





Solid, Polycrystalline Sources 

In Chapter 3 we discussed a theory due to Abragam and Pound 
(14535) related to the correlation spoilage in polycrystalline 
sources. We will now apply this theory; the justification for 
the application will be discussed in Chapter 6. 

We proceed by obtaining the observed anisotropy, corrected 
for positrons but not for geometry, of the correlation as a func- 
tion of (letw where (Lo is the angular frequency associated with 
the maximum electric quadrupole or magnetic dipole level split- 
ting, i.e. the energy difference between the highest and lowest 
levels in the previously degenerate state which has been pertur- 
ted by the interaction. From the plots of the above functions 
(Figs. 5.4.4-2,3) we can obtain a lower limit for.ts for each 
assumed interaction. Then knowing an upper limit on ta, the 
nuclear lifetime of the intermediate state, we can set a lower 
limit on flo, 

Let A''!' be the observed anisotropy corrected for syste= 

matic errors as defined in Section 5.35.2, 

A'' ‘be the observed anisotropy corrected as above but 
also corrected for positrons, and 

A be the observed anisotropy corrected as above but 
also corrected for geometry. 


The anisotropy is defined: 


(1) >. * Ameo) 4 
wO90) 


If we express the correlation function in a LeGendre poly- 


nomial expansion: 


(2) wie) = 1 + AodPo(cos®) + A,P,(cose) 


then the anisotropy is given by: 


a 





Laie 4-" om % 
8 — 4A5 + Shy 

From (AI-36): 
(4) Alt = attr_g 

where G is the positron correction term. 

In Section 5.2.3 is was shown that for sufficiently small 
values of € and J, the geometrical corrections cause: 

(a) no mixing of the coefficients for finite solid angle 
effect, and 
(b} only negligible mixing of the coefficients for finite 
source effect. 

We can then write: 
(5 ) Wit! = 1+ FohePolcos@) + Fydgh4(cose) 

wnere the F's are simple reduction factors arising from 
the above mentioned geometrical corrections. 

Then by (5) and (3), the observed unspoiled anisotropy 
corrected for positrons becomes: 
ee A! lel te Naa 

8 — 4F odo t SF yA, 
And from the spoilage theory mentioned above, the spoiled 


enisotropy as a function of Qots, defined above, is given by: 


(7) A''( O.%0) = 12PoGo(N. te dbo + SF4Gq( etn dAg 
8 — 4FoGo(a.tw Ao - F464 (notw Jay 
where the G's are expressed in Chapter 3 as a function of 
Wety where Wo is the angular frequency associated with the mi- 
nimum level splitting in each case, and F,4,are the undisturbed 


coefficients corrected for positrons but not geometry. 
p 


2 


We restate the expression developed in Chapter 3 fo 


er = 





G's for the two tyves of interactions in volycrystalline sources 
assuming classically describable fields of axial symmetry. 
(a) Slectric Quadrupole Interaction; 
(8) Go (twetw) = 4372 , _.057 4 +545 __ 4 +228 - 
T+ Gd" Ter 1+ QW. TH 1416 we Cw 


G wots Y= | Feo Pbead 0 254 ~ 095 


+ 


-{- 7 ee ee 
14 Godin 149wset  1416Wetw” 


(>) Magnetic Dipole Interaction: 


ey G_ (wwetv) = (2 4s 2A a 1 
a L+ Wer tr 7 1 4 theses 


Ga, (wet) = 111 + £222 Ly as 
1+ We ty 1+ 4ué Tw 


a i 
+ Tyoun 7 TF | 


where the subscripts "e" and '-" refer to e-ectric and 
magnetic interactions respectively. 

We now wish to express‘Llo in terms of Wo , both of which 
nave been defined above: 
Electric Quadrupole Gase: 

It can bexshown that the minfmum level splitting in an 
axially ratte vec field for integral I is given bys 

a-V 
ealSe)3 

fh Cale 


where e@ is the quadrupole moment as defined in Chapter 3 


(10) Cape — 


2 
and (2 | is the gradient of the electric field along the axis 
of symmetry. 
The maximum splitting is given by: 
(Se) 
_ Cait) SE 
(11) Loe Ee 
h 4(2T-1) 


ey ee 





So that the ratio: 
(12) (lee = Woe rs (For integral I) 
Magnetic Dipole Case: 

The minimum level splitting here is just the angular ve- 


locity associated with Larmour precession: 
ja. H 
, t 


(13) Wom = 
where #4 is the magnetic moment, and the maximum split- 
mone 1s$ 
(14) SlLem = 2H (For integral or half 
f integral 1) 

So that the ratio: 

m5) (Lom — Woe 21 (For integral or half 
integral I) 

We use Slcemwas an argument instead of wombecause as the 
fields take on lower degrees of symmetry than axial the angular 
frequency associated with the maximum splitting,-tlom, to a first 
approximation remains constant, while this is not necessarily 
so for Wom. 

ror illustrative purposes only we also include in our 


plots (Figs. 5.4.4-2,3) the previously discussed (Chapter 3) 


square wave G' tg, 





k 
aK S 
* a oe eeees K 
(16) Gi (ate) == Le Wetw e2ginh( Gian) 
20 
e Wein — [ 


where for I=2 we have (again from Chapter 3): 
(17) 5S Sr. +9484 = 4 

where the S's are the half widths in radians of the 
square waveSe 


The purpose of including (17) and (18) is to show graphically 





now G's of a widely different nature than those based on condi- 
tions ef axtal symmetry will yield curves of approximately the 


same shave near the origin (where Notw—> 0) 


Liquid Sources: 

The analysis of the spoilage of the anisotropy in liquid 
sources as a function of the viscosity of the source shows that 
we cannot apply Abragam and Pound's (1453) liquid theory (Chap- 
ter 3). The reason is as follows. Abragam and Pound develop 


coefficient attenuation factors for liquid sources in the form: 


| 
( + \k CN 


where Nk is directly proportional to the "correlation time", 


f—is) Gk = 





ee , of the liquid. Correlation time can be thought of as 


the average time it take for an ion in a liquid to interchange 
its neighbors. As Cc gets much larger than @n we approach 
the condition of the solid state; hence for the liquic theory 
to apply tc.<tn. We observed however, no detectible dimunition 
in the anisotropy of the liquid sources until Ce was of the 
order of 3x107” seconds which is at least three orders of mag~ 
nitude larger than Cn 3; hence we cannot apply the liquid 
theory to this spoilage. This point will be discussed further 
in Chapter 6. We present therefore the liquid data mereiy in 


tabular form. 





Part 5.4 - Results 
Section 5.4.1 - General 
In this section we will consider only the results of the 
experiment and ieave their interpretation to Chapter 6. 
Lifetime measurements were taken on the intermediate state 
of qqail4 by the method of delayed coincidences using a coinci- 
dence circuit similar to that of Bell et al. (1852). 
Retult: Qwe< 2:3 IO a S.. 
Section 5.4.2 - Gorrections 
Decay Scheme Corrections 
Calculation of 6 (AT=31) 
Following are entries for equation (AI-31) for the 
calculation of the positron correction, 6 « Only sige 


nificant plus and minuses are shown. 


a — 2 we fin; = tees Joay = 1632 
[Lg = 0.65 Pesg = 0649 Jisg = 0657 
Sq 227916) Sop = 5014168) S_ = 2.9616) 
S500 50(5) $54, = 40205) S, = 5.84(5) 
Pea a2. 4246) S55 = 501916) S, = 2.78(6) 
SOg (90) == 1.27(3) Note: Numbers in paren- 
Sip (180) = 2.79{3) theses indicate 10 to 
Ciep(90) =m 3,04(2) the power. 
Con (180) = 6.19(4) 
a, ==0.039 + .046 “ Xp == 4(-5) & (+5) 
Eye = noes + Foca eee) 


Ws mena 2,48(-2) 
* Second approximation did not change 


Entering these figures in (AI-351) we obtain: 


0 — 0.015 * .004 


Dis 





Geometrical Corrections 
Finite Source Correction 


For significance see Appendix [1I. 


€. = .010 ¥. = .007 


Then (ATIIe7) becomes: 





Co = 0.972B5 + 0.013B5 + 0.009B, 
Co = 0.942B5 -++ 0.0843, 
Cy 0.893B, 
Finite Solid Angle Correction 
Se@:Bo. (5.3.5-13) and Fie. 504.2-1: 
Qo = 1.1900 Q4 = 1.6885 


Section 5.4.3 - Undisturbed Correlation 

General procedure 

a. Performed least squares fit on points corrected for sys- 
tematic and decay errors as defined in Section 5.35.2. 

be Applied positron correction to points and obtained ano=- 
ther least squares fit. 

ec. Corrected the LeGendre coefficients for finite source 
effect. 

d. Corrected the LeGendre coefficients for finite solid 
angie. 

e. Introduced the closely observed and corrected "“anisoe 
tropy oniy™ measurements in comparison with “d” and obtained 
final coefficients by using center of "cverlap” of the two 
areas in Fig. 5.4.3-1. 

These results are tabulated below. 

wie) = 1+ apcos%e + agcost@ = 14 AoPoleose) 


+ A,P,lccse) 
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Procedure ao a4 Ao Ay 





a .039 +.040 .087 1,041 

b .059+,040 .052+.039 .067+.035 .012 .009 
¢ .o70+.037 .013.010 
rl .0447-,030 .1007.075 .0sst.043 .o2et.oie 
e .055+.062 .099t.070 .o90+.o11 .o2et.o16 





Experimental points together with curves "b" and “e" are 
shown in Fig. 5.4.32. 

We list below the possible correlation coefficients for 
various spin assignments in the cos“~9 scheme from (1H40) for 
pure multipole transitions of order L. The spin assignments of 
the three states involved are given without parentheses; the 
multipole order of the transition is given in parentheses. 
Since cqil4 is an even-even nucleus we assume the ground state 
has a spin-parity assignment of zero-pvlus and the first exci-e 
ted state two-plus; hence for the second transition the only 
possibility is electric quadrupole, E2. 


Group a ~ First Transition Dinvole - Second Transition 


“Quadrupole: 7 
| =e A 
1(1)2(2)0 -0.353 O 
2(1)2(2)0 0.429 O 
8(1)2(2)0 -0.103 O 


Group b = Both Transitions Quadrupole: 


0(2)2(2)0 3000 4,000 
1(2)2(2)0 5.000 -5.3350 4 
2(2)2(2)0 -1.154 A eetl 
3(2)2(2)0 O -0.353 
4(2)2(2)0 0.125 04> 


75 = 
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As will be shown in the table which follows, none of the 
observed coefficients supports any possible spin assignment in 
Group ae In Group b the 4(2)2(2)0 agrees with Steffen (1551) 
but none of the others: however it is closer to the others than 
any other assignment in Group b. However, Steffen and Zobel 
(1852) have shown subsequent to Steffen's first measurement of 


this correlatiom (1951), that the second excited level (1278 


114 


KEV) in Ga is fed directly from the ground level (one-plus) 


4 


of Int* by Kecapture (See Fig. 4-1) which indicates a spin 


of two for this level in ceee The assignment of one-plus to 


114 35 based partly on the election rules 


for the beta decay to the ground state of aa. This proceeds 


the ground level in In 


by an allowed transition (log ft=4.5), and since the ground 
state of oa“ has a spin assignment of zero-plus (being an even- 
even nucleus), the ground state of a must have a spin-parity 
of one-plus or zero=-plus on the basis of the selection rules. 


4 
The K-capture transition to the second excited state of catt ’ 


being allowed, implies that the spin of this state (in cqtt*) 

is zero-plus or one-plus or twoeplus. If it were zero-plus, 

the first transition would have to be E2 and the cqail4 cascade 
would be 0(2)2(2)0, with no possibility of modification by multi- 
pole mixtures. And as will be seen, this fits no observed cor- 
relations hence, we can eliminate zero-plus. Of the two remai- 
ning possibilities, one-plus or two-plus, the former is ruled 
out because it 1s not possible to find a multipole mixing ratio 
which will fit any of the observed correlations, whereas as will 
be seen, this is possible with a spin-parity assignment of two- 
plus, yielding the following cascade: 2(1,2)2(2)0. In addition, 
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in assigning two-plus to the second excited state of Cd » we 


cee 





are supported by the knowledge that most other even-even nucleii 
have either two-plus or four-plus for this level. 

We find that we can obtain a good theoretical fit to the 
observed curve by considering the first transition to be a mix- 
ture of Ml and E2. From the tables of Ling and Falkoff (1L49) 
we construct Fig. 5.4.5—3 which gives the LeGendre coefficients 
of the correlation as functions of the quadrupole mixing ratio. 
In Fig. 5.4.5.-3 the curve Ao is part of a tilted eccentric el- 
lipse and A, is a straight line; only that part of Ap correspon- 
ding to a relative phase of zero between the quadrupole and die 
pole matrix elements is given, for the best fit is in this req=- 
gion. Using Fig. 5.4.d=3 we obtain the mixing ratio and resul- 
ting theoretical coefficients. These are given below in compari- 
sor with those of Steffen and Zobel (1852) and Johns et al. 
(1354): 


ge (%) M1 (%) 
Steffen and Zobel 4.4 95.6 
Johns et al. 4.0 96.0 
Daubin 4,2+.006 '9538-b. 006 


Following.is tabulated the coefficients as observed by 


various workers: 


Name Ref. Ao a, Ay A, 
Ae Klema & McGowan 1K52 084 ~ 106 elll 025 
iD: Johns et al. 1J54 -~,0028 e190 «087 e059 
Co Steffen 1951 ~ 125 e042 SOS ~ 009 


d. Steffen & Zobel 1852 090 «070 . 096 O15 
eC. Daubin © 055 ,coo » O90 A OE 


+ eee! + "Forno = Joni + one 
ee Theoretical 086 064 »090 014 
2(M1,E2)2(2)0 
E2 —_ 4.2 % 
Ml = 95.8 % 


With the exception of "c" these curves are shown in Fig. 5.4.3=<4. 
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Section 5.4.4 = Disturbed Correlations 
Polycrystalline Sources 
See Section 5.3.7 for a discussion of how the following curves 
are constructed and the meansing of the different paramters used. 
The source of InCls dissolved in a dilute aqueous solutiom 
of HCl yielded the true, or undisturbed, correlation. This will 
be discussed further in Chapter 6. Since we use the anisotropy 
so measured as a starting point for the curves in Figs. 5.4.4-2,3, 
we list below the apparently undisturbed anisotropy together with 
the correlation coefficients: 
A't(0) = .124 T .009 
AS'(O) = .074 F .037 
BO) = (Ollsias geome 
The observed anisotropies for the various polycrystalline 


sources are given below: 


Source Descriptiom Ai! 
(1) InClz Salt 097 +,013 
(2) Frozen (~-104°9C) dilute © Uo ae (Ona! 


aqueous solution of In0l. 
in HC1 (,4N) 
(3) Frozen (a/-104°C) dilute 5087 1s B10 
solution of In0l. and HCl 
in > 99% glycerin, 
(4) Average of (1), (2), (3) 2.095 1 007 
(5) Average of (2), (3) 095 +,007 


In Fig. 5.4.41 we show the theoretically predicted per- 
turbed anisotropy based om Abragam and Pound's theory of clas- 
Ssically describable, axially symmetric perturbing fields and 


an electric quadrupole interaction (1453) (See also Section 
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5.3.7 and Chapter 3). The purpose of this plot is to show the 
Overall shape of the function. The predicted anisotropy A!', 
corrected for positrons but not geometry, and normalized to the 
value actually obtained in the undisturbed case, is plotted a- 
Bainst the magnitude of the interaction,tlotw, where_Nlois the 
angular frequency associated with the total quadrupole splitting 
(Eq. 5.3.7-11). Also shown is the "hard core" anisotropy to be 
expected from fields of rhombic symmetry for I=>=2. In this case 
the “hard core" G's are Go= G4 =2/7 (1A53) 

In Figs. 5.4.4-2,5 we again plot theoretically predicted 
anisotropy normalized to the undisturbed value versus magnitude 
of the interaction wsing the same total splitting angular fre- 
quency discussed above. In both Figs. 5.4.4-2 and 5.4.4—5 we 
assume axially symmetric, classically describable fields; in 
the former we assume an electric quadrupole interaction and in 
the latter a magnetic dipole interaction. The elliptic looking 
curves drawn within the statistical limits of the various ob- 
served points represent approximate lines of equal probability 
for one standard deviation, i.e. 68% of samples should lie within 
closed curvess. The lower vertical tangents to these curves re- 
present the lower statistical limits of the interaction,Notw . 
The observed anisotropies of sources (1) and (5) as defined 


above are considered. 
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Liquid Sources 
We tabulate the results of the observations om liquid 


sources. See also Fig. 5.2.2-2. 


Temperature Viscosity Mean Correlation Alt At} 
(Deg. Cent.) (Centipoise) Time (Seconds) 
207k 3 zi (2m) .139 | .124 
+ .008 + .009 
20% 3 800 8(-9) Ass © wae 
~ ewe +. 010 te..008 
Bagel «dk 4000 4(-8) 0133 6118 
+ 800 +> 4010 tc Jema 
~14t 2 33000 3,5(-7) .107 , ,092 
= 10000 + ‘008 + +009 
~28t 2 600000 6(-6) 121 .106 
++ 200000 + ico a em. 


Note: Numbers in parentheses under "Correlation Time" ree 


present 10 to the power indicated. 
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Chapter 6 = Interpretation and Discussion 
Part 6.1 - Undisturbed Correlations 

One of the questions frequently asked when a value of ani- 
sotropy is alleged to be from the “undisturbed correlation" is, 
"How do you know the correlation to be undisturbed?" We shall 
attempt to answer this question using two lines of argument. 
The first involves a comparison with the results obtained from 


ita ial 


Cd » and the second involves a method of successive approxi- 


mations among the results pertaining to cqilt, In both arguments 
we shall assume that in the dilute aqueous solution Abragam and 
Pound's theory of spoilage in liquids applies; this is to say 


that we will assume that @c< ‘Cw although we do not know this 


to be true with certainty in the case of cqilt, Both arguments 


make use of (3-26) which give the spoilage G, factors in liquids 
as a function of the electric quadrupole moment, lifetime, and 


correlation time. 
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In the comparison of cqtll with ca in the aqueous solu 


tion, we begin by quoting from (2H53) a value of 500 mes for 


“wes 
Cb 32) where eqiil is the quadrupole moment of the first 


nls of 





excited state-of Cd « We make the further assumption that 

2 
i | is the same in the oa??+ aqueous solution as it is in the 
gqil+ aqueous solution. We find that when these values are put 


. -10 
into (3-26) with the additional value of Tw S oe CeeO + seconds 
for the Cal!* 1ifetime and an estimation of 10721 seconds for 
the value of Cc, that we have but negligible attenuation, 


1.e. Go ~G, ~~ 1. Further, attenuation in the aqueous solu- 


4 
tion does not become appreciable in cait unless we make the 


unrealistic assumption that we have a quadrupole moment in Ca 


iL 
about thirty times as large as that found in eqtt - We sum 
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marize by saying that on the basis of this test, our aqueous 
solution gives the true undisturbed correlation if we assume 
that the ratio of the quadrupole moments of call4 to caltt is 
no larger than thirty to one. 

To begin our second argument let us assume that out re- 
sults from the dilute aqueous solution at room temperature repre- 
sent the true unperturbed correlation. Then we use the aniso- 
tropy so determined as a starting point for calculating the va- 
lues of the interaction in the polycrystalline wear GS), 
5.4.4-2,3). Then, let us take the magnitude of Ct \32* 
determined and place it in (3+26), just as we did in Mt first 


111 interaction. If this in 


argument with the results of the Cd 
turn causes a reduction in the G's from unity then we will ap- 
ply this correction to the original anisotropy and re-calculate 
the interaction. If, however, such a substitution leaves the 
G's essentially equal to unity, then we can say that our aqueous 
solution does indeed represent the undisturbed correlation. We 
have also made the implicit assumption, which seems justified 
at least in order of magnitude, that the average electric field 
gradient experienced by the nucleus in liquid sources is no lar- 
ger than that experienced in polycrystalline sources. If we 
make these substitutions we find again that the G's again re- 
main essentially equal to unity, and from this we conclude that 
we may consider for our purposes that the aqueous solution yields 
the truly unperturbed correlation. 

In regard to the complete correlation, it is interesting 
to not that our value for the mixing ratio in the first transi- 
tion: M1 —> Yess = .o06 

E2 —> 042 + .006 
agrees within the limits of error to those ratios found both by 
= o6 — 








Johns et al. (1754) and Steffen and Zobel (1852) (Section 5.4.3) 
and thus lends added confirmation to the accuracy of these values. 
On the other hand our fit to the theoretical curve based 
on the above mixing ratios (Fig. 5.4.3-4) seems better than that 
of the others. This adds credence to our values for the LeGen- 


Gre coefficients of the correlation: 





Ao = “090 2 .o11 


A .o22 + 016 
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Part 6.2 - Disturbed Correlations 
Section 6.2.1 - Polycrystalline Sources 

Our results for the magnitudes of the interactions and mo- 
ments are based on various general and particular assumptions. 
We shall first state the assumptions, then our conclusions based 
on these assumptions; then we shall discuss both the validity 
of the assumptions and the firmness of the conclusions. 

General Assumptions: 
1. The perturbing fields are classically describable. 
2; The perturbing fields are static. 
Se The perturbing fields are axially symmetric. 

Particular Assumptions: 
a-l. The interaction is electric quadrupole in nature. 
a-2, The interaction is magnetic dipole in nature. 
bel. The primary perturbing field is due to the crystalline 

structure and the average electric field gradient at the 


x6 statvolts/cm2. 


nucleus is about 1.6x10 
be? The primary perturbing field is due to the electronic 

cloud of the atom of which the nucleus is a part and the 

average electric field gradient at the nucleus is about 


435x101 statvolts/om“. 
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bed The primary perturbing field is due to the crystalline 
structure and the average magnetic field at the nucleus 
is about 10* Oersteds. 

b~4 The primary perturbing field is due to the electronic 
cloud of the atom of which the nucleus is a part and 
the average magnetic field at the nucleus is about 


10° Oersteds. 





Zs s&s 5 
Conclusions: C) eA _ ie 
: vo = & z AY 
et AvVe \, 
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Let Adm = gpl a ° = ee 


where eQ is the nuclear quadrupole moment, 





His the nuclear magneton, and 
g¢ is the nuclear gyromagnetic ratio. 
Then we list below for the various sources, a lower limit 
for the AV and the moments based on our known upper limit for 


10 


the nuclear lifetime ( < 2.3x10° seconds), and the statistical 


limits from Figs. 5.4.4-2,3: 


Source: Partic. Assump. 


(a) Magnitude of the Interactions: 


(1) Incl, Salt a-1 Atoe >, 700 mes 
(2) Frozen Solutions of InCls a-1 Athe yy 990 mes 
(3) InCls Salt a2 AVom > 163 mes 
(4) Frozen Solutions of InCl.z a=2 AVoom 2 214 mes 


ae 





(b) Magnitude of the Moments 


(6) InCls Salt a-1,b=-1 la | »y 269% 
(6) InClz Salt a-1,b-2 |Q | 77 ee) 
(7) Frozen Solution of In€lz a-1,b-1 lal >; «8s 
(8) Frozen Solution of InClz a-1,b-2 aul Yass 
(9) InClz Salt a-2,b-3 | pl % 42 
(10) InCigz Salt a-2,b-4 ) pl 7 -42 
(11) Frozen Solution of InClz a-2,b-3 ne >, 56 
(12) Frozen Solution of InClz a-2,b-4 | jel o56 
Note: The units for the quadrupole moment: x107"* em”. 

The units for the magnetic dipole moment: Nuclear Magnetons. 
Discussion 


That the fields at the nucleus are classically describable 
would hold in both the magnetic and electric cases if we could 
neglect the excited condition of the cadmium atomic configuration 
following K capture; for the ground atomic state of cadmium is 
an § state. But from the discussion in Chapter 3 we see that 
we cannot negiect the “excited configuration effect" whose fe- 
neral theoretical interpretation requires the introduction of 
the quantum mechanical coupling between the nucleus and elec- 
tronic configuration. This in turn would be complicated by the 
uncertainty of the total electronic angular momentum, J, fol- 
lowing K capture. However, as pointed out by Abragam and Pound, 
in the case of the electric quadrupole interaction (particular 
assumption a-1) we can still use the classically described 
fields because in polycrystalline sources orbital electronic 
This is not possible in the 


degeneracies have been removed. 
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case of the magnetic interaction because of the remaining spin 
degeneracy. This effect can be seen on semi-classical grounds 
as follows. When a nucleus is acted upon by an external field 
it also reacts upon the field; i.e. just as it is possible for 
an external field to induce nuclear transitions, it is also pos- 
sible for the field produced by the nucleus in a transition to 
induce transitions in the configurations giving rise to the 
initial field. If the energy levels of the states giving rise 
to the initial are well separated (no degeneracies) then the 
nuclear transition reaction can cause no transitions among 
these states; this corresponds to the classically describable 
field. However, on this basis it is easy to see that the 
reaction due to a nuclear transition could induce transitions 
among degenerate states which are the sources of the initial 
field, These “reaction induced” transitions among degenerate 
states could, in turn, alter the initial field. In such a situ- 
ation the only valid treatment is a quantum mechanical one. 
With this limitation clearly in mind, we have made certain cal= 
culations (assumption b=-4) in which we neglected this quantum 
mechanical cotpling in the case of magnetic interactions with 
the electron cloud, realizing that there still remains a pose 
sible spin degeneracy even in polycrystalline sources. However, 
in an order of magnitude estimate, as this calculation is, this 
procedure seems justified. In the light of the above, conclu~ 
sions (9)=(12) should be viewed with a certain amount of reser- 
vation. 

That the perturbing fields are static seems to be well ta- 
ken, (1H52) indicates a shell relaxation time of greater 


than 107° seconds for ionic crystal impurity centers which is 


age. = 





some two hundred times as large as our nuclear lifetime. Hence 
the catl* nucleus experiences essentially a constant perturbing 
field due to the "quasi-static configuration" effect (Chapter 3). 
That the perturbing fields are axially symmetric is not at 
all clear. As discussed in Chapter 35 even if an individual 
crystal did provide axially symmetric fields (where two of the 
field gradients are equal) at the nucleus these fields would 
probably be distorted by inter-crystalline pressures in a poly-= 
crystalline source and by nuclear recoils out of axially symme- 
tric sites following K capture. However, a redeeming feature 
is that the hard core values for rhombic fields (where all three 
field gradients differ) are not much different from the hard 


cores for axially symmetric fields; as we saw in Chapter 3 for 


wa. s 
Axially Symmetric Rhombie 
G,(1im) 13/755 10/35 
G,(1im) 29/63 16/63 


And from Fig. 5.4.4-1 we see that the hard core anisotropies 


in the two cases are? 


7 A(lim) 
Axially Symmetric © 046 
Rhombic »035 


The Gy, ( {Lot ) for axially symmetric fields fall off quadratically 
with iL.tw asSiitw —> 0, with a second derivative nearly equal to 
the square of the angular frequency associated with the total 
splitting (Le) times the nuclear lifetime ( tw ) (3-21). 
Simple arguments indicate that this same general quadratic bee 
havior with the same relation to total splitting should apply 


for fields of lower symmetry. Since the total splitting is 
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approximately independent of the exact degree of field asym- 
metry for a given value of maximum field gradient, we would ex- 
pect that the G for axially symmetric and the G, for rhombic 
fields would lie very close together for small values of the 
argument, (lotw , 

In regard to particular assumptions a-1 and a-2 we will make 
no effort to favor one type of interaction. Abragam and Pound 
(1A53) show that even if a magnetic decoupling experiment shows 
negative results, i.e. no detectable increase in anisotropy, 
as reported by Albers-Schonberg et al. (2A53), this is not con- 
clusive to the absence of a magnetic interaction. We will thus 
keep a-l and a-2 on an equal footing and not try to argue for 
either. However, if the interaction is electric quadrupole 
(a-1) our quantitative results are in general more credible than 
otherwise because of the discussion above relating to classical- 
ly describable fields. 

As for particular assumptions bel to be4 relative to the 
sizes of the fields; these are only order of magnitude approxi- 
mations and are not to be taken too literally. This indeed is 
one of the gréatest obstacles to the more fruitful use of dis-= 
turbed angular correlations for the determination of moments, 
i.e. our ignorance concerning the magnitude of fields experienced 
by the nucleus. The value relating to the electric field gra= 
dients in crystals (b-1) came from the appendices of (1G53) by 


15 


averaging a number of cases which ranged from 8.5x10 TO 5e5X 


ic statvolts/cm”. These values of the gradients were cale 


10 
culated from instances where both Q and the splitting in megae 
cycles were known. The value relating to the electric field 


gradient due to the electron cloud (b-2) is from (1T49) and 
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is based on the gradient due to the p-electron of iodine; this 
would correspond roughly to one of the first excited states in 
cadmium. The values for the magnetic fields (b-3,4) are from 

(1H52) and (1F52). 

The only conclusions backed by direct experimental evi- 
dence relative to spoilages are (1) = (4) relating to the magni- 
tudes of the splitting frequencies; these are to be taken as 
firm within the limits of statistical validity and of the indi- 
cated assumptions. 

Conclusions (5) - (12) concerning the moments themselves 
are uncertain due to the obscurity as to the sizes of the per- 
turbing fields as mentioned above. However, it seems we can e- 
liminate conclusions (9) and (11) as being far too large based 
on typical nuclear magnetic moments; then we can say that if 
the interaction is magnetic it arises from the "quasi-static 
configuration" effect, i.e. an interaction between the nucleus 
and the magnetic moment due to its own electron cloud. For 
comparison purposes the quadrupole moment of the first excited 


state of pal0® 


» an even-even nucleus, as estimated by Abragam 
and Pound (1453) is at least .2x10~<* om, 

The complexities of the processes going on in the source 
introduce a great deal of uncertainty in the results, and it 
must be borne in mind that an exact theoretical calculation of 
the "quasi-static configuration" effect may well lower our 
values of the interaction. However, if the highly excited state 
of the electronic configuration following K capture is primarily 
responsible for the perturbing interaction, as it seems to be 


in the case of Steffen (19853), then these large fields could 


provide an excellent means for the determination of the moments 
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of very short lived ( < 107° seconds) excited states. There 
seem to be three great experimental hurdles at present in the 
way of specification of these moments more precisely: 

(a) Measurements of nuclear lifetimes less than nor’ secs. 

(b) Measurements of the magnitudes of the fields experienced 

by nucleii under various conditions. 

(c) Specification of the nature of the interaction, i.e. 

whether electric or magnetic in nature, 

Furthermore, supposing we could perform (a) and (c), but 
not necessarily (b), then comparative measurements of the spoi- 
lage due to the “quasi-static effect" would give us the ratio 
of moments for different isotopes of the same element. Or, 
using known moments, the spoilage ratio in two isotopes could 
give us information as to the atomic physics of the "quasi- 


static" field, 


Section 6.2.2 — Liquid Sources 

The data from the variable viscosity glycerin solutions 
yields no information as to the magnitudes of the interactions. 
This is because (as can be seen from the table in Section 
5.4.4) at the viscosity at which spoilage may perhaps be ta- 
king place (about 3x10° centipoises) the theory of spoilage 
in liquids (Chapter 3) developed by Abragam and Pound (1453) 
no longer applies. At this viscosity the liquid correlation 
time is about 10°. seconds which is some 700 times as large 
as the nuclear lifetime, and a basic postulate of the liquid 
theory is that the correlation time be shorter than the nue 
Clear lifetime. At such a large correlation time, our nucleus 


with a lifetime of less than or equal to 2.3x1071° seconds, ex- 


oe 





periences practically static fields very similar to those asso- 
ciated with the solid state. 

Another phenomenon in the liquid sources that bears inves- 
tigation is that there is no observable spoilage in the range 
where (OT wS‘CeST00%; here the nucleus is in a region of transi- 
tion where it is begining to experience fields which are prac- 
tically constant, and hence related to polycrystalline pheno~ 
menae A possible answer may be that in liquids even of this 
high a viscosity there is some mechanism which shortens the 
atomic shell relaxation time appreciably, and thus reduces the 


average field at the nucleus. 


Part 6.3 =- Statistical Validity 

In arriving at the statistical validity of the results we 
will make the following tests: 

1. A Chi-Square test on the average anisotropy of the 
four sources, viz., aqueous solution at room temperature, In0l.z 
salt, frozen aqueous solution, and frozen glycerin solution, 
all considered simultaneously. 

2. A Chi~Square test on the average anisotropy of the 
last three of the above mentioned sources considered together. 

3S. The two-point statistical test developed in Section 
5.2564 wherein, if one assumes two points come from the same 
normal distribution, the probability of finding them separa- 
ted by more than so many standard deviations is given. When 
worked in reverse, i.e. given two points not closer than so 
many standard deviations, it gives an upper limit on the pro- 
bability that they are from the same normal distribution. This 
gives only an upper limit due to the possibility of their 


being from different distributions. We will perform this 


=a G4 w 





test, comparing the aqueous solution with each of the three 


other sources. 


Results of the Tests: 
Test No. Probability that a random sample gives no 
better fit: 
1 041 
2 5S, 
Upper limit on probability that two given 


values are from same distribution 


) Aqueous Sol. - InClz Salt - 0352 
Aqueous Sol. ~- Frozen Aq. Sol. - - 089 
Aqueous Sol. - Frozen Glyc. Sol. - 001 


None of the above tests is in itself conclusive. The fact 
that Test 1 gives a probability of .041 is certainly not con- 
clusive evidence that the aqueous solution and the other three 
sources do not give results which really lie on the same distri- 
bution. However, from Test 2 one can state with a high degree 
of certainty that the three polycrystalline sources give re- 
sults which lie on the same distribution or on distributions 
which are volte: 1 Oke together. When, on the other hand, one 
considers Test 1 in the light of Test 2 and the fact that a 
spoilage correlation actually seemed to exist, i.e. any de- 
viation in the anisotropy from the aqueous solution was in 
the downward direction and this under conditions that predict 
"spoilage", then one has strong reason to believe the exis-— 


tence of a true effect. 
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Appendix I - Positron Annihilation Component of Coincidence Rate 

At 180 degrees the Inti4 cascade gamma coincidence rate ig 
‘poliluted" with a rate due to the positron annihilation component 
of the decay. We wish to determine what fraction of the observed 
true coincidence rate at 180 degrees is due to this annihilation 
component. 

We will consider three different sources, which we will de- 
Signate as follows: 


Nee Inv4 


extended source 
Ni - Na“© extended source 
N, = Man point source 

Notations: 
G) ~- effective counter solid angle 
E€ e- counter efficiency 
- channel factor, i.e., fraction of detected 
pulses falling within "window". 

The probability that a single gamma ray of energy E, emitted 
isotropically, will register a single channel count in the ith 
channel is given bys: 

NEC = Wee ee Sei 

Note that all three factors are energy dependent. 

A more precise definition of (w& and hence of € will be 
presented later (App. II); for present purposes, however, we 
may speak as though the counter had a uniform efficiency € for 
all gamma gays within solid angle \W , and zero efficiency oute 
side ° 

We designate the Inl14 branching ratios as follows: 

Xp = positron emission 
Xk = electron capture for the process which 


leads to the 722-556 KEV cascade in call, 


= Ges 


The intensities of the various gamma rays are as follows: 


ofx Ne - Each of the two tnii4 


cascade gammas 
ZL Xp No ~ Both of the Inil4 annihilation gammas 
Ni - The Na<® cascade gamma ( Cs cascade) U=1)2 


oe annihilation gammas (2152 


2Nu -=- Both of the Na 
In this discussion we will consider gamma rays of four dif- 
ferent energies. For simplicity in notation each energy will be 
represented by an integer as follows: 
l= 722 KEV 
2=~- 556 KEV 
5 - 511 KEV 
+ 


1277 KEV 


Single Channel Rates 


114 be Soi 


Det the Gth (lL =G.b ) channel rate for In 
Let the & th (¢:ayb) channel rate for Na“< be Sib J=ly2 
Thens: _ 

(1) Sot = No | Xx (Mii +2) + Zoleac| 
but sinceXp<< Ak, (1) becomes:: 

My S0i & Ak Ne Chie + h2i) 


and similarly for Nae”: 


(3) Sac = Ni Chae trac) 8 eh? 

In order to simplify the following derivations, we will use 
a term we call “tne” single channel rate, Se. In order to pre= 
sent the relationship of Se to Sea and Sob , consider the case 
where there is present only a gamma ray of energy J . Then by (2)' 
me Soi = Xe Na Nic i gir 

If we define So: 


(5) Se = 4 aac Sais 


= OF «= 





and Nie 
(6) ee, = Ok No Ai 


Then, by (5) and (6): 
(7) \i = lie Nie 
Extending the above to the case where there are two gamma 


) 
rays, we have a more general definition for the 8 $ 


(8)) So = Ax No GN a 


Coincidence Rates 


es savers gamma coincidence rate be Cog (0) 


Let the In 
Let the Inti4 annihilation gamma coincidence rate be Cop (180) 

Let the "a annihilation gamma coincidence rate be Cip(ix%) s=!y2 
Let the Na@* cacade-annihilation gamma coincidence rate be Cigp (80) 
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The fn cascade gamma coincidence rate may be expressed as 


follows: 


(9) Cos (0) — Xx No W(e) | Madzh + Kahle} 


Combining (2), (5) and (8), we see that:: 


(10) Che #2 \ = (\ia + Nee | 


If we introduce (7): 
Ga) 2\.e = Viren + do dIb 
Then we reals’ 
(12) Cog(@\ = 2X« Nea We) bike 
The annihilation gamma coincidence rate at 180 degrees, 


however, is expressed: 


\2 ab 
(13) Cop ({20) = _2%p No £2 haa hab 
O33 
where (lL is a reduction factor due to the finite size of 


the source, i.e., (L=1 if the mate of every annihilation gamma 
entering one counter, enters the other counter. One obtains this 


with a centered point source, but not with an extended one. 


- GB = 





By (7). we have directly: 
2% No Ne 


(14) Cop(iso) = 
== 


Since out knowledge of Nais limited at best, we can be rid 
of this factor by considering the ratio of (14) and (12) at 180 


degrees: <e 
C180) eet Na 
Cog, C180) ok WUBI, NENZ 


We are rid of Noagain, and have another relation, if we 
consider the ratio of the single channel to coincidence rates 
at 90 degrees where we know there is no positron component. Then 


by (12) and (8): yd 
(16) So me i 


Cog eGo. = _ 2 \N (90) wie 


We now define the term: 


(17) Meas = <= 
And for the _ of completeness: 
; = War : 5 
(18) L321 = PS i cll =: Che 
2 Mu 


Since one can actually measure [Lis 2\ by comparing relative 
areas under lines in the single channel spectrum, it is worth- 
while to show the explicit relation between (18) and (17). 

From (11) we have: 
es ay, | Sepok — 

- >} ae 
And introducing (7) and the definitions (17) and (18): 
= 
(20) fpa= lo| aya + ps2 | 
or, the arithmetic mean of the single channel pes ° 
Then, by (17) and (16): 
Se 


\ a, —— | + 
ae Cag (90) 2 W (30) AZ 
By substituting (21) into (15), again using (17), we come to: 
(29) Cop(i80) __ 4 ol p W (40) CQ. So \3 
Coq (180) Xk W Lid) (4 Cog (so) ( I+ 20) 


If we compare the single channel to coincidence rates for the 


Na@* sources, both point and extended, we find an expression for; 
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5 wal @ Meum 152 NSC) S2. 
C 2e (180) Si 


We get a grip om \3 by considering the coincidence rate 
in Na<® between the annihilation and cascade gammas, which for- 
tunately is isotropic. By a procedure similar to the above we 


obtain: 
(24) ha = 14 SEE (14 2 poe] 


By substituting (24) and (23) into (22) we have: “ 
(25) Copliso) _ , op W (9a Cie citer ig p go) Sa Sa (1+ A peat) fa 


TY NM Lo a a renee 
Cog (igo) H&K W (140) Coa Crp (120) Cog (40) Si (1+ [har] 


All terms have been expressly defined except CU3 which is 
defined more precisely in App. II, where also a method for mea- 
suring it is derived. 

wR 

The presence of the second power of W (Jo) in (25) is unfor- 

tunate, but not seriously so. It would be much neater to have 
W 90) = 

the ratio Yi gpPPear To this end, we will proceed to cal 

culate W (40) ; 

(26) W (a) — =n (l+aArcos*@ + CL up cost @) 

where KX is a normalizing factor whose magnitude is deter- 
mined by the following constraint: 
(27) Wieldus = 417 

Unit SPHERE 

If we carry out the Le PO tees it turns out that: 
(28) Kk = [| { /ap ean eas a4] 

"ai" and "qQ," are to be determined by a least squares fit 
to the observed data. 

By (28): ot, 

| 
(29) W(9o)} = K = [1+ nz + 15 a4] 
And substituting (29) into (25) we have: 


ae re AL 
(30) Capliso) ale Wark CyplisdCiap (Io) Jo Sa ((4 2 34 | Ma 
Coalited 4x Wiig) (t33 Cop (Ko) CO S6.(90) 57 mas bray i 
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Let us define and O as follews: 


A K Cip (0G 1gp (40) Bo S82 phat (\+ Mae)” 
jis 


4X% Cox Cap Ciney Gaae, (qo) SP C (+ Ma 


Con C10) 


Then (30) becomes: 


(33) fr = Wao) 
W (180) 


And since:: 


mes) C.o(@) = Con (0) +-Coe(@) 


where Cal@)is the total, true rate, it is easily shown that:: 


(35) Wcigod  __ see ( | 


Wao) = Cotso) \I+ 


And by (32): 


WW (1\<0) <= Col1go) _ Q 
W (90\ -_ Colo) 











(36) 


We thus arrive at the pleasant result that the correction for 
positrons to the observed anisotropy consists merely of subdtracting 
some number which we can calculate explicitly by (31). The pro- 
cedure must be carried out by successive approximations, however, 


because of the uncertainty in KK in (31). 





Appendix II ~ Effective Solid Angle 

In equation (AI-25) we find the quantity ©) 3 , which has 
been referred to roughly as the “effective relative solid angle”. 
It is the purpose of this discussion to obtain for this term a 
precise functional form. 


As a precise definition we use an equation similar to AI~14,. 


2 
my Cay (B01 = 2N2 hs 
Ws 
(a 
| 
(2) Wa — 2N2 NZ 
C2P(w) 


But we can represent x as the integral of some relative 
response fine tong GM, over the counter, so normalized that: 


~ = ae F (x) das 
(3) 47 
> eS 
where OQ is the angle subtended at the source between the 
detector axis and the element of solid angle, AW ; 


We can write the coincidence rate at 180 degrees for the 


No (Na, point) source as follows: 


pa 
me Coe utc) = ae \F ty ED 
: a aa (6 ae 


Substituting (4) into (2) we haves remembering (3): 


a | ( Fide 
" | § F “va dew | 


We observe an angular resolution curve, (7 (6), from the di-~ 


(5) (ie 


rectional correlation of the annihilation radiation. Making use 
of the counter axial symmetry of the function F (x) and refer- 


ring to the diagram on the top of the next page, we can express 
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Counter 





(+(@) as follows: 
(6) G(@) = « | Fie) F(t) dao 


(7) 


where X=O(@,6,)\ and K is an undetermined constant. 


We obtain the resolution curve as a relative function, 


Ne) mere gy NS Fey Fades 


H(@) = = 
(>(0% \ F(A) du 


We can expand F (0) in LeGendre Polynomials:: 


FQ) = mm Cm ee (cos ©) 
v 





and using the addition theorem: 


Faget = Cm | Pe (cos 0) Pm(cos®) 
rep ae fib Pe (605.0) Panccs VCs 04] 
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The terms in cosud vanishes on integration over @ . Hence: 
+ 


(10) [ \ Fe) dus| H (@) _ am \ Fonts Pm (x) Paamddx >; (Cos @) 


Making use of the orthogonality of Pn(x): 


- - 2s 
5 a | eae 1s 
(11) | ai Cerdeo| H(@) = 27 2. ( = Pm (cos @) 


Multiplying through by Pn(cos@) and integrating again we 


come down to: 


ee fl 2 
(12) ae = a | | \ Hie) Pn teoss eos. X 
From (8) we have: " k ( F*ce)du| 


+1 
4s) +1Co = \F(@) des 


Or, using (12): 


4 [ ( Fredo] = an | ( Woerd caso 6 Fede | 


Then by (5): +] 


(15) WwW, = = ) Herd cose) 


Notes A factor increasing the apparent source size for 
annihilation Padiation is the range of the positrons. Although 
this is taken into consideration by (15) it may be of interest 
to know the positron “halo” around the source. In the decay of 


In-t4 


the positrons have a maximum kinetic energy of about 1 
mev. This means that their range in the lucite source holder 
is about .13 inches, as compared to length 1.10 inches and dia- 


meter .o7 inches which geometrically define the extended source. 
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Appendix III - Geometrical Corrections For A Finite Source 

When one observes a directional correlation function, one 
is observing a phenomenon in the macroscopic domain. The ob 
served function consists of a weighted sum of micro-scopic func- 
tions, each related to an element of source and an element of 
solid angle on each counter. The function, W(6), which theory 
predicts is related to the microscopic domain, i.e. it is con- 
cerned with well defined directions in the decay of a single 
nucleus. However, the observed function, which for the pure 
poses of this discussion we will designate W(®),results from 
radiations proceeding in directions uncertain to the extent of 
the counter solid angle and which result from the decays of mil-~ 
lions of nucleii whose position is uncertain to the extent of 
the source. 

Let our observed function be: 
(1) W@) = 2_01Pi(cos®) 

where (# ve the angle between the axes of the two counters. 
Note that we do not restrict the i's to even integers. 

And let our well known micro-scopic function be: 
(2) wie) =. 1+ ApPo(cose) + AgP4(cose) 

where 9 , as usual, is the angle between the propagation 
vectors of the two radiations involved in the cascade. 

We now ask, "What is the functional dependence of the 
C's on the A's?" 

To answer this we must first specify the model of our de~ 
tection system very carefully. For our first model refer to 


the diagram on the next page. 
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Detection 


ag 






Sphere 






om 
Geometrical 


Extent of 
Counters 


X Counter Resolution Curve, € (ot) 


In the diagram dn, represents an element of source, di), and 
dus, are elements of solid angle on the counters. The counter 
angular resolution curve for one counter is shown, the radial 
direction being the coordinate of relative response. 


We make now a few comments about this model. Certain 
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artificialities have already been introduced in that we are re= 
presenting our counters by spherical caps concentric with the 
Origin. Obviously, this is only an approximation, since a real 
counter has q radial extent and all detection does not take place 
at exactly the same distance from the origin. We assume in this 
model that the relative counter sensitivity, E(od) , is symme- 
tric about the counter axis and is independent of the slight de- 
viations from normality to the counter surface of radiations from 
different parts of the source. This, it seems, is a good approxi- 
mation for small sources, i.e. where the direction of the propa- 
gation vector from various parts of the source to an element of 
counter surface, does not change appreciably in direction. With 
these comments in mind we can then write our observed correlation 


function as follows: 
(3) W(@) = W (Eid Er (ud dme dw. dw 


Wee 
4% An SourRcE < eC a 
Now, in principle it is possible to integrate (3) and thus 


obtain the functional dependence of the C's on the A's (1,2). 
We would then be making our solid angle and finite source cor= 
rections simultaneously and accurately to the limit of validity 
of this model. But in practice this integration is exceedingly 
tedious. 

We now proceed to alter our model slightly to facilitate 
the integration (3). We do this by separating the corrections 
for finite counter solid angle and finite extent of source. 4 
As was pointed out in Section 5.3.5 the corrections for finite 
solid angle with a centered point source cause no “mixing” of 
the coefficients. In this case, assume that our observed func= 


tion is W'(0), where: 
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(4) wi(e) = 1 + BoPol(cos®) + ByP,4(cose) 
Then it can be shown that: 


(5) Boy = Roxon 
where Ao, is the undisturbed coefficient and Qo, is a simple 


reduction factor, i.e. no mixing of coefficients occur. 


If we now alter our model and assume that each element of 
source can be first corrected in turn for solid angle without 
introducing mixings we can separate the two corrections, i.e. 
for finite solid angle and for finite source. This is equivalent 
to saying that in our model the spherical caps can be tilted 
slightly to bring each element of source in turn to the cap's 
spherical center, without deviating too far from the true state 
of affairs, 


Then from this model we can construct another diagram: 


Oo dome \ 
= 12> i. 


e: . 











ene. 


—_ 
i =? Counter #2 
ys — 


ie — 


Counter #1 





Here we have already carried out the solid angle integration 
so each counter can be replaced by a spherical point at its cen- 


ter: 
(6) W on 


expresses the observed correlation. If now we assume Ri = Ro, 
and that the source is in the form of an upright right circular 
cylinder centered at the origin, of “half-height" z, and radius 
ro then we obtain an explicit expression for the C's (1) in 
terms of the B's (4). Letting z9/R = ae , and o/h = S » we 


have after carrying out the integration in (6) to second order 


in C. and iy 


Zu ce 
F ) Co sao By “te E (Bo + Be +84) AG (289 ~ Bp — 84) 


— g(a “3. — 5B) + bt Bot B,) 
| i el g, ( 3,) + a - 3Bp— 384) 
. a Ge | “Fo we Je 3) 
= oe —3B,) + 3 —2B,) 
i ~ 5B, ) 


We can express the B's in terms of the A's by (5). We thus 
have answered the original question, iJe. we have expressed the 
O's in terms of the A's. Our answer is accurate to the limit 


of validity of the second model discussed above. 
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